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ABSTRACT 


This thesis demonstrates a method to determine optimum, tactical movement paths 
for a specified vehicle and or small umit based on the operator’s cognitive decision 
processes, as well as the physical effects of terrain and environment on mobility. The 
approach uses psychometric techniques inherent to the Generalized Value Svstem (GVS) 
to determine a “Power Function” based on a specific tactical scenario and given equip- 
ment configuration and provides a means to determine the Tactical Movement Potential 
(TMP) for each terrain cell. This cognitive value in an interval scale can then be trans- 
lated into the same scale as the physical continuum using techniques proposed bv L.L. 
Thurstone and W.S. Torgerson. The cognitive time value based on the user’s decision 
process 1s then added to the phvsical traversal umes for each cell computed from output 
provided bv The Condensed Army Mobility Munagement Svstem (CAIMS). This ren- 
ders a value mapping which can be optnnized using one of several existing algorithins. 
The Dikstra Algorithin is used in this demonstration model. The resulting sets of path 
points are optimized for speed time and the cognitive tactical considerations evaluated 
using these psychometric methods. Lhe movement path and resulting times can be used 
in combat planning and modelling. This output ts also particularly nmportant in deter- 
nuning the tume values needed to compute the Situationally Inherent Power (SIP) of the 
GVS. This methodology could be apphed to almost any tactical decision process in the 


development of expert svstems and models. 


THESIS DISCLAIMER 


The views expressed in this thesis are those of the author and do not reflect the of- 
ficial policy or position of the Department of Defense or the U.S. Government. 

The reader is further cautioned that certain vehicle system input parameters used 
and portions of the computer program developed in this research are not valid for all 
scenarios of interest. While every effort has been made, within the time available, to 
verifv all computer programs: they cannot be considered validated. Anv application of 


the methodology using these programs 1s done so at the user’s own risk. 
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I. INTRODUCTION 


A. GENERAL 

The ability to accurately portray the battlefield environment is a key requirement in 
combat model development. This same ability 1s required in the construction of tactical 
decision aids and tools. One major factor to be considered in any of these applications 
is the effect of terrain, vegetation, and weather on the mobility trafficabilitv of weapons 
systems and units. Another major factor is the effect of these same parameters on Visi- 
bility, detection. and acquisition by various weapons systems. While the preceding fac- 
tors are objective for the most part and can be accurately measured; other factors are 
not so clearly defined and may possess large variances. This is certainly the case where 
a human decision process 1s required such as the decision logic used by various combat 
unit and vehicle commanders in tactical movement. A clear distinction must be made 
between mobility as modelled in relation to the physical constraints of the battlefield and 
movement which requires interaction with the user decision maker or must be modelled 
taking the pertinent cognitive considerations as well as the physical constraints of the 
battlefield into account. 

In any case, all of these factors and more must be accounted for in our combat 
models and tactical decision tools. There are a number of existing high resolution 
combat models and tactical decision aids which incorporate many of these functions. 
Most of these packages use one of several available versions of digitized terrain data and 
either actual, historical, or user input weather data. Most combat models also incorpo- 
rate appropriate tactical decision logic algorithms to control some key functions such 
as detection, acquisition, and other essential combat actions of player units. Some tac- 
tical decision tools also do this. However, at the present time there is no combat model 
or tactical decision tool which combines the appropriate decision logic to determine the 
Optimum path a unit or vehicle should traverse in order to simultaneously maximize 
speed and trafficabilitvy, minimize inter-visibility and detection, and appropriately bal- 
ance their effects in terms of the decision logic. 

The need and desire of commanders and combat planners to have and use such 
tactical decision tools and to incorporate them into combat models is of definite benefit. 
Some commanders presently have a tactical decision aid in the form of the Army Mo- 


bility Model (AMM) and the Condensed Army Mobility Management Svstem 


(CAMMS) which provide trafficabilitv. speed. and trme contour mappings for use in 
route selection and other tactical decisions. There are also any number of combat 
models available for study where a commander or combat analyst can select appropriate 
terrain data bases and manually input desired unit or vehicle path points and or posi- 
tions. This is the case with present JANUS based models, the Battalion Combat Out- 
come Model (BCOM), and others. However, there is no tactical decision tool or combat 
model that can systemically determine sets of optimum, tactical path points for the 
commander or for use in analysis. This is also true of almost all models where tactical 
decision processes are represented. 

The first requirement in developing such a model is to properly identify those vari- 
ables which could andor would affect the tactical movement or other decision process 
such as the mission, threat. corresponding locations, ranges to Various threats, vehicular 
agility or localized speed, allowable completion time, and others. [tis soon obvious that 
a large number of variables, many having probable interactions with other variables, 
exist. [The possible combinations or instances is quite large and each variable may vary 
significantly between instances. Structuring this maze in order to be able to extract usa- 
ble data on the cognitive process is no small task. 

‘Another requirement is finding a theoretically sound method for evaluating large 
groups of variables simultaneously and then being able to compute a specific value 
Which appropriately measures their combined effect on the process. This field of study 
has a name all its own. Psvchophysics. 

A third requirement ts to find a theoretically sound method that allows for the 
mapping of these psychophysical values into an appropriate physical continuum so that 
they mav be used concurrently with phvsical attributes and variables in the model. Thus 
ficld of study also has a name all its own, Psychometry. 

The art of combining and applying these techniques is known as Psychometrics. 
We will use these techniques in conjunction with an existing mobility model, CAMMS 
(Version 2.0), to provide total movement time values based on physical engineering 
traversal times and the translated psychometric times for each terrain cell and subse- 
quent traverse segment on the battlefield. These values will then be used as network cost 
values in a form of Dikstra’s Algorithm which renders the set of minimum cost optimum 
path points for that scenario. 

The organization, purpose, scope, and assumptions of this research are explained in 
this Introduction. Background information on previous route selection methodologies, 


terrain representations, physical attributes of terrain considered in mobility, cognitive 


considerations Which impact on movement, and the desire to place such technologies in 
the hands of the commander in the field is presented in Chapter II]. The development 
of the “Psychometric Function” using psychometric procedures inherent to the Gener- 
alized Value Svstem, GVS, is shown in Chapter III. Application of these techniques and 
construction of the model algorithm 1s detailed in Chapter IV. Results of the algorithm 
unplementation on an AT-386 personal computer and limited verification of these results 
are presented in Chapter V. Proposed enhancements to the existing model and areas re- 
quiring additional study are presented in Chapter VI. Finaily, Chapter VII discusses 
possible utilizations for this methodology and proposed applications with respect to 
tactical decision and training tools, research and analysis tools, combat models. and the 
Generalized Value Svstem. 

It is of significant interest to note that prominent Soviet military analysts, Nikita 
Silomseyev in particular, have proposed a modei framework based on decision theoretic 
developments using such state variables and decision functions since the 1970's. Two 
Specilic areas discussed have been tactical decision processes and movement of the line 


@mecontact. (Ref. 1: pp. 21-26] 


B. PROBLEM STATEMENT 
!. Purpose 
This thesis examines and demonstrates a methodology for determining opti- 
mum, tactical movement paths for a specified vehicle and’or small unit based on the 
operator's cognitive decision processes, as Well as the physical effects of terrain and en- 
vironment on mobility. The resulting sets of path points are optimized for speed time 
and the cognitive tactical considerations evaluated using psychometric methods thereby 
rendering a truly optimum. tactical movement path and optimum traversal time. The 
movement path and resulting times are to be used in combat planning and modelling. 
This output is also particularly important in determing the time values needed to com- 
ailtestne oittiationally Inherent Power (SIP) of the GVS. 
2. Scope 
This research will be limited to one specific scenario of interest in order to 
demonstrate the methodology without undue involvement in various scenarios with dif- 
ferent sets of fixed state variables. The specific scenario used is discussed in detail in 
Chapter III. This research will not involve detailed studies of relative effective range used 


in the algorithm and discussed in Chapter IV. There will also not be any extensive ver- 


ification and validation of the results of the computer application of the algorithm due 
to a lack of time and resources. 

This research will examine the use of a movement model based on the method- 
ology for a single vehicle, unit of specified type and mission against a threat array of 
specified tvpes performing a certain mission while in a selected area of operations and 
known terrain. 

3. Assumptions 


|. The mathematical forms of range dependency equations studied by Bonder and 
Farrell and used to portray perceived lethality as relative eflective range Cumecme- 
lated to those found in JAY Soare cCnrect: 


2. The survev results, interval scales. and regression analysis are valid tor themmanrae 
of values surveved only. Hlowever, these should encompass the vast majority of 
Instances: 

3. The translated cognitive time scale and traversal time scales are additive andmmen 


multiplicative. 


4. This methodology ts applicable to all types of tactical decision processes when the 
appropriate state variables can be identified. 


Il. BACKGROUND 


A. PREVIOUS ROUTE SELECTION AND TACTICAL MOVEMENT MODELS 
[. The Dynamic Tactical Simulation (DYNTACS) Model 

The DYNTACS Model was developed by the Systems Research Group at Ohio 
State University from 1964 to 1969. The model is extremely high resolution in order to 
predict individual unit performance in armored combat engagements from both the de- 
sign and operational perspectives [Ref. 2]. Dynamic programming is used to run and 
solve the simulation which consists of a driver program and 34 subroutines. At least 15 
of these algorithms are related to or affect tactical movement in the model [Ref. 3]. An 
explicit surface terrain representation of the battlefield is used in order to determine 
plane departure points for a predetermined grid size [Ref. 3: pp. 57-66]. This represen- 
tation is closely related to and can be derived from exphcit grid terrain data available 
from a number of sources [Ref. 4: pp. 3.3-3.15]. 

The model determines an optimum route for an advancing element or unit by 
optinuzation where “Tactical Difficulty” is mininuzed. Tactical difficulty is computed 
Meme a heuristic of the form, /D =(1 + E)T , where E is the difficulty computed due to 
a set of factors along each route segment and T is the travel time for that segment [Ref. 
3: pp. 94-116]. Such heuristic procedures are computationally efficient but do not guar- 
antee an optimal solution. The set of factors and corresponding values for the difficulty, 
E, were obtained using comparative judgement techniques. Travel times are computed 
using engineering models for each weapon system. It 1s interesting to note that difficulty 
Peaeiunetion Of the effective range, actual range, and threat disposition [Ref. 3: pp. 
97-100] and that psychometric methods similar to many we will examine were used to 
obtain these values. The major shortcoming with the DYNTACS approach 1s that the 
cognitive and physical scales are related multiplicatively using a heuristic due to existing 
computer linutations rather than developing a truly optimal technique. 

2. The Simulation of Tactical Alternative Responses (STAR) Model 

The STAR Model is similar to DYNTACS in most respects except that dynamic 
programming is discarded and the optimization is performed in FORTRAN applving 
Dijkstra’s Algorithm. The additional concept of the “Sliding Pattern” is used to restrict 
the search for the path [Ref. 5: pp. 23-34]. The route selection is performed sequentially 


from start point, to a horizon or horizons, to an objective by sliding the optimization 


along the network. Unfortunately, the identical heuristic equation found in DYNTACS 
is also applied in this model [Ref. 5: p. 33] rendering it suspect. 
3. Other Route Selection Methodologies 

There are any number of similar movement models proposed over the past 25 
vears and certainly all of them cannot be mentioned. However. it is necessary to intro- 
duce another type of model which is probabilistic in nature. One such model was inves- 
tigated and formulated by Sitmourang [Ref. 6] in 1981. This dynamic route selection 
model performs an optimization in order to minimize the probability of being attrited 
enroute to the destination. The attrition probabilities along each route segment are de- 
rived from two sources. These are enemy elements capable of destroying the tactical el- 
ement and terrain factors which end movement. Speed reduces the probability of 
attrition due to the enemy but increases the probability of attrition due to terrain. espe- 
cially rough terrain [Ref. 6: pp. 21-22]. The total probability of being attrited. P. is taken 
to be the sum of the probability of being killed by the enemy, PKFOE, and the proba- 
bilitv of being killed due to terrain’speed, PKRVEL, which is minimized using Dikstra’s 
Algorithm [Ref. 6: p. 34]. 

There are a number of problems with this methodology. The model assumes 
perfect information in computing probability of kill due to enemy. This is certainly not 
the case in most tactical decisions. The probability of attrition due to terrain speed 1s 
computed using a heuristic equation with a power factor or constant for all combina- 
tions of conditions which is applied to a ratio of recommended speed to actual speed. 
This value is then used multiplicatively with normal expectancy in order to determine a 
probability of attrition due to speed in given terrain [Ref. 6: pp. 24-26]. No reference or 
theoretical support for this technique is provided. This certainly makes such a method- 


ology suspect. 


B. CURRENT MOBILITY MODELS 

Mobility models provide comprehensive mappings by terrain cell or other construct 
of the measure or ability of weapons systems to physically traverse types of terrain in 
varying conditions. Mobility models should not be confused with movement models 
which attempt to describe how tactical elements traverse the battlefield based on both 
physical and cognitive considerations. The two primary mobility models used today are 
The NATO Research Mobility Model (NRMM) and The Condensed Army Mobility 
Management Svstem (CAM MS). Most mobility routines used in models today are based 
on this family of models. The NRMM and CAMMS are both derived from the Army 


puooiity \iodel (A S11) developed in the 1970's. The CAMMIS is currently used as the 
base inodel for: 

1. The Planning Analvsis Work Station (PAWS) 

2. The Terrain Analvsis Work Station (TAWS) 

3. The Airland Battlefield Environment (ALBE) 

eee ine Digital Terrain Support System (DTSS) 


The model is also used by numerous Army, Defense, and other government activities to 
include The Defense Mapping Agency, The U.S. Army Command and General Staff 
College, The United States Military Academy, The Army Research Institute, The U.S. 
Army Engineer School, The U.S. Army Training and Doctrine Command (TRADOC), 
and The TRADOC Analysis Command (TRAC). The CAMMIS is used as the base mo- 
bility model in this research due to its applicability and acceptance in the modelling 
commuuuty.[Ref. 7] 
I. The Condensed Army Mobility Nlanagement System (CAMIMS) 

The CAMMS vehicle performance model was originally developed by simplify- 
ing the existing AMM through addition of a vehicle preprocessor and restricting the 
model to one vehicle type at a time. Obstacle performance is also determined statistically 
based on vears of testing and data collection. This allows the CAMMS to operate on 
personal computers versus larger mainframe computers for which the AMM was devel- 
oped. The model is much too large and detailed to be discussed fully here; however, an 
examination of the physical aspects considered in the CAMMS, output available, and 
current uses of the model is worthwhile. 

a. Physical Considerations Affecting Mobility 

The CAMMS off road prediction model evaluates three general categories 
of terrain data. These are surface geometry data, surface composition data, and vege- 
tation data. Surface geometry data includes slope information, surface roughness factors, 
and all obstacle data. Surface composition data primarily relates to the tvpe of soil and 
its strength which are greatly affected by weather. The vegetation data includes stem 
diameter, tvpe, and spacing data along with visibility information. These data are used 
in a number of routines to compute surface traction, resistance, obstacle traction, inter- 
ference, gap crossing capability, vegetation effect, and ride dynamics which are combined 


to produce a highly accurate speed prediction for the vehicle.[Ref. 7] 


b. CANIMS Output Mappings 
The CAMMS 1s widely used throughout the Army and particularly in the 
European Theater of Operations. A range of mappings or Tactical Decision Aids. 
TDA‘s, are currently available in the CAMMS. These include: 
I. Vehicle Predictions 
a. Speed 
b. Time 
c. Reasonslemnted 


d. Comparison of Two Vehicles 


2. Foot Soldier Prediction 
3, Maneuver Damage Prediction 


4. Soil Strength Prediction 


5. Route Analvsis Prediction 


The vehicle predictions and especially speed and time mappings for the digitized terrain 
grid are of most significance to this research [Ref. 7]. These will be addressed concerning 
their potential use in subsequent chapters. 
c. Technology to the Field 

The CAMMS takes advantage of recent advancements in the personal 
computer field to provide an extremely powerful yet portable and inexpensive tool for 
the commander in the field. The most current CAMMS requires only an IB\{ compat- 
ible AT personal computer with a math coprocessor, DOS 3.2, EGA graphics, | mega- 
byte of RAM, and sufficient hard disk drive storage in order to operate. All them iyms 
mentioned previously are available in a matter of minutes provided the appropriate 
digitized terrain data bases are available. The only shortcoming in the current system 1s 
the lack of standarization in terrain data bases and lack of data in general for some re- 
gions. The CAMMS is being used extensively as an aid on REFORGER exercises and 
as an operational planning aid in some tactical units. 

The methodology developed during this research will also be derived So as 
to make it conformable for use on a personal computer. This will allow the model to be 


used as a portable tactical decision tool in the hands of the commander. 


C. COGNITIVE CONSIDERATIONS IMPACTING ON MOVEMENT 
The discussion up to this point makes it quite apparent that tactical movement 


modelling is much more involved than mobility modelling. The major question becomes, 


mow do we model the cognitive decision process? There is no definitive answer to this 
question and numerous techniques are available and feasible depending on the applica- 
Meme It 1S apparent that the mission, terrain effects, and threat all influence the 
combatant or tactical decision maker to some degree. The decision concerning where and 
when to move 1s made based on the combatant’s evaluation of how the enemy can affect 
his mission accomplishment and ‘or survival as he moves to each possible location. There 
is also a question of the trade-off made by most humans between surviving and accom- 
plishing a mission within a certain time. Sonie people are more risk averse than others. 

The process 1s synonomous to shopping for fruit or vegetables of different types and 
deciding which items to purchase based on their size, shape, color, texture, and price. 
The shopper makes a determination when comparing different items as to which 1s best 
in terms of overall quality relative to all others and what he will purchase with a certain 
amount of funds just as a combatant makes a determination of the potential lethahty 
facing him based on the types of enemy, their range, and their capabilities based on 
range and disposition in order to decide what risk he can accept when moving on the 
battlefield. 

An enumeration of the variables identified as relative to tlus tactical movement de- 
Sienprocess and how they are structured for use in this research is provided in detail 


ie hapter {{I. 


Ii. PSYCHOMETRICS AND PSYCHOPHYSICAL ANALYSIS IN THE 
GVS 


A. GENER SIE 

The Generahized Value System, GVS, is an axiomatic value system used as a tool for 
evaluating the power or potential of entities at present and future times. The formal 
structure of GVS was first proposed by Professor Arthur Schoenstadt of the Naval 
Postgraduate School in 1985 [Ref. 8: p. 3]. This procedure outlines a method of deter- 
mining the state of an entity over time. allows us to compare the values of entities rela- 
tive to one another, and is consistent for all entities on the battlefield. The wales 
assigned to each entity are referred to as “power” in previous GVS documentation. 

Power is a subjective value given to each element on the battlefield and combines 
both inherent and derived power. Derived power consists of the power an element pos- 
sesses due to its ability to influence support other friendly elements [Ref. 9: p. 33]. In- 
herent power is the ability of an element to directly affect the outcome of a battle [Ref. 
9: p. 32]. Inherent power may be categorized™as Basic Imherent Power sien 
Situational Inherent Power (SIP) which are the power an element has at full strength 
and the power an element is predicted to have at some future time, respectively. The 


SIP is highly dependent on the available time. 7, , which is determined by the timesdo:= 


main networks forthe respective battleieids (cls mace 

The methodology proposed in this research directly supports use of the GVS by 
providing a means of determining actual optimum time values for the terrain grid net- 
work representing the battlefield. The methodology will also provide actual optimum 
time values and sets of path points for use in other combat models using conventional 
Lanchester attrition or firepower scores as long as they utilize explicit terrain grid net- 
works for representing the battlefield. It will also provide excellent tactical movement 
planning and decision tools when incorporated into the appropriate mobility models. 

Modelling the cognitive process in order to determine appropriate values for a given 
tactical movement scenario and instance and to then be able to map these values into 
the physical continuum, time, on the battlefield requires an understanding of how atti- 


tudes are measured, how to develop respective interval scales, how to convert cognitive 


values in interval scale into physical values in a ratio scale, and how to solve for the 


“Psychometric Power Function” which allows the user to solve for this value given anv 


set of realistic variables at different levels. 


B. THE MEASUREMENT OF VALUES 

The correet measurement of psychological values 1s fundamental to evaluating in- 
stances or combinations of variables and subsequently being able to use them in any 
physical construct. [lowever, there are numerous wavs to develop measurements and 
Seemeresults in different scales with Varving characteristics. The three major character- 
istics of measurement are order, distance, and origin. Distinctions between tvpes of 
measurement and types of scale depends on what combinations of these characteristics 
est ief. 10: p. 15}. 

A nominal scale exhibits none of these characteristics and will therefore not be 
evaluated in terms of potential use [Ref. 10: p. 17]. An ordinal scale does possess order 
and may have an origin; however, the scale does not exhibit distance which results in 
Mo ability to discriminate between instances in terms of a value [Ref. 10: p. 16]. This 
limits potential tvpes of usable scales in psychophysical analysis to interval scales. 

I. Interval Scales 

Interval scales alwavs possess the properties of order and distance. It 1s also 
often possible to determine a rational origin for such scales [Ref. 11: p. 196]. Interval 
scales Which exhibit all three major characteristics are known as ratio scales. These scales 
are the most important scales used in math and science and the term “measurement 1s 
Siem pestricted to these tvpes of scales {Ref. 10: p. 31]. Interval scales can be determined 
in a variety of wavs for sets of instances and scenarios. Three powerful methods are 
Continuous Response Scale Judgements, Comparative Judgements, and Categorical 
Judgements. 

a. Continuous Response Scale Judgements 

A continuous response scale allows a judge to indicate the value of an in- 
stance on a continuous scale between a minimum and maximum bound, normally 0 to 
100. This method lends itself to Analvsis of Variance (ANOVA) testing commonly used 
to determine relationships. The scale value of each instance is easily obtained and the 
summary statistics along with the distribution of responses for each instance are also 
easily obtained. However, one problem with this approach is that judges have difficulty 
actually determining values to assign between various instances [Ref. 8: p. LO]. This 1s 
magnified the larger the number of instances becomes. For example. only five variables 


broken down into five distinct levels would result in 5° or 3125 unique combinations of 


I] 


variables requiring responses. It 1s extremely difficult for a judge to realistically place a 
value on each with this number of instances and this approach was rejected due to this 
factor. 
b. Comparative Judgements 

The Law of Comparative Judgements presents the judge with all possible 
pairs of instances in order to determine which has the greater value in each pair. This 
discriminal process relies on the principle of “just noticeable differences” developed by 
L.L. Thurstone. The results of these judgements can be tabulated to provide a discrimi- 
nal difference and deviation [Ref. 11: pp. 39-47]. Application of this theory is done using 
the method of paired comparisons. A matrix solution based on the frequencies of the 
responses 1s used to compute an interval scale. However, with n instances it requires 
oe comparisons [Ref. 10: pp. 166-179]. For the five variables broken down into 
five levels, this results in 3125 instances and 4,881,250 separate comparisons. This is far 
more comparisons than 1s feasible in most studies and too labor intensive in terms of 
data collection. 

c. Categorical Judgements 

The Law of Categorical Judgements parallels the Law of Comparative 
Judgements in many respects. The judges evaluate each instance independently and se- 
lect a rating category they feel best represents the general value of that instance. Ihe 
categories are understood to be a mutually exclusive set of successive intervals which 
collectively exhaust all possible values of the continuum. Each rating category has de- 
scriptors and explanations which serve to help the judge with his task [Ref. 12: p. 1]. 
Normally distributed responses for each instance across the sampled population 1s a key 
assumption of this method. This is easily observed by examining the frequency matrix 
in Appendix B. Computing the interval and or ratio scale is relatively straight forward 
using a matrix solution [Ref. 10: pp. 205-239]. A clear, concise method 1s given by Pro- 
fessor Glenn Lindsay. Discussions with Drs. Parry and Lindsav at The U.S. Naval 
Postgraduate School assured that such an approach would be theoretically sound and 
vield usable results without requiring excessive survey and data collection effort and time 


once the scenarios Were properly structured. 


C. THE EXPERIMENTAL DESIGN 
As previously addressed in the Problem Statement, the purpose of this study is to 
develop and demonstrate a methodology for determining an optimum, tactical move- 


ment path based on both the cognitive and physical considerations or constraints of the 


battlefield. As was also stated in the Introduction, devising and structuring tactical sce- 
narios Which insure the correct combinations of state variables or instances impacting 
on movement logic are addressed to the appropriate sample population can be an enor- 
mous and complex task. Administering the data collection effort once scenarios are de- 
veloped is another major step in the experiment. 

Wie process Occoimes more structured once the data are obtained®The Law of 
Categorical Judgements is used in this case and a Well defined method for obtaining and 
translating the interval seale from the responses is available. Multiple regression analvsis 
can then be applied to determine the functional relationship of the scale values for each 
Instance to the set of state variables. The resulting regression equation becomes the 
Mes cnometric Power Function’ needed to determine values for use in GVS and other 
models. 

I. Scenario Development 

Identification of all possible state variables or those items needed to describe the 
specific tactical situation is the first step. This was accomplished in several ways by ex- 
amining previous and ongoing tactical movement studies, reviewing doctrinal informa- 
tion On movement. and through discussions with groups of selected officers and combat 
modellers familiar with tactical movement. The following 14 variables are identified as 


the significant factors affecting decisions on tactical movement: 


I. Mission 


eee lime available 

3. Equipment and resources available 

4. Threat equipment and capability 

5. Threat intention or nussion 

6. Range to the threat 

fee OVEr 

oaee Olicealment 

9. Environment, both weather and obscurrants 


10. Area of operations or theater 

11. Speed or vehicular agility 

12. Range or distance to the objective 

13. Obstacles, both natural and man-made 


14. Artillery 
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The problem now becomes how to organize these variables into scenarios con- 
taining organized combinations or instances in order to provide a means of measuring 
certain attitudes concerning tactical movement. The number of state variables and their 
respective levels must be structured to reduce the required number of responses while 
maintaining the desired degree of fidelity in the experiment. 

One technique ts to eliminate any variables not significant to a majority of 
sources. [his was the case with the variable artillery since the major impact doctrinally 
concerned obscurrants which fall under environment. In addition, other sources consid- 
ered artillery to influence movement only through the use of improved or special 
munitions which usually fall under obstacles. 

A second technique is to fix selected state variables for a given scenario, espe- 
ciallv those which are categorical variables and not suited to a regression analysis. The 
respondent or judge is made aware of these in the definition of the scenario and the re- 
sulting scale and function are only valid for that specified combination of fixed variables. 


The following state variables are fixed when using this methodology. 


Le Siission 


tJ 


Threat intention or mission 


Area of operations or theater 


din TOS) 


Equipment and resources available 


in 


Threat equipment and capability 


A third technique 1s to identify which variables may be combined when actually 
considered by the decision maker. This is the case when evaluating cover, concealment, 
environment. and range. All of these combine to determine the degree of inter-visibility 
or “line-of-sight” between the element and the threat. This line-of-sight or LOS and the 
number of lines-of-sight are what become critical to an element manuevering on the 
battlefield. There are also a number of excellent line of sight algorithms available for use 
in combat modelling. 

This same technique is applicable in the case of localized speed or vehicular 
agilitv. Obstacles, environment, and vegetation combine with the terrain itself to deter- 
mune speed. An element moving over the terrain is not concerned with soil moisture or 
tree stem spacing but with the total effect all these variables have on his speed. The 
combatant evaluates movement based on speed because he knows that the greater his 


speed, the less likely it is he may be engaged and the less time he will require to reach 


14 


his objective in general. The Condensed Army Mobility Management System, CAM MS, 
provides excellent and highly accurate speed mappings for digitized terrain data. 

The variable “range to threat” also requires further investigation. A combatant 
must evaluate many threats simultaneously and they are often different tvpes with 
varying capabilities. The combatant makes his tactical movement decision based upon 
the greatest known or suspected enemy capability which can influence his nussion at that 
moment. This becomes his greatest perceived threat. As mentioned previously, the 
combatant makes a determination of the potential lethality facing him based on the 
types of enemy, their range, and their capabilities based on range and disposition. A 
variation of the range dependency equations used in many Lanchester models will be 
feeceto calculate the “Relative Effective Range’ or RER of each enemy system and al- 
low the analvst to compare the lethahty of various systems at different ranges much the 
Same as the combatant. These have several mathematical forms studied by Seth Bonder 
and Bob Farrell from 1969 through 1974 [Ref. 13: No page]. The nunimum RER of all 
threat systems with line-of-sight at a particular location then becomes the range which 
equates to that range given in the survey and used in the model. 

One last observation concerns the variable time. Time actually becomes the de- 
pendent variable for which we wish to determine a relationship. The range to the objec- 
fives also directly related to and affects the time domain network inherent to GVS. 
These two variables are not included in the set of state variables given for each instance. 

This now reduces the number of state variables to the five fixed variables iden- 
tified previousty and the folttowing three key variables. 


Peleoculized speed 


tJ 


Number of lines of sight 


eyeeicanece to the greatest threat 


The three kev variables must now be broken down into distinct levels in order 
to construct the set of instances needed for the opinion survey. An extract of the actual 
survey used 1s located at Appendix A. A conscious decision was made to keep the re- 
quired number of responses at a level which allowed the the survey to be completed in 
an hour or less by a respondent. The variables are each divided into four levels which 
cover the majority of conceivable occurrances for each variable. This results in 4 or 64 


instances to be evaluated and which must be appropriately organized for the survey. 


2. Survey Construction and Execution 

This is only a pilot survev given to a limited population of U.S. Army and Ma- 
rine Corps officers assrgned to the U.S. Naval Postgraduate School, TRAC- Monterey, 
and Fort Ord. A total of 62 surveys were distributed during a one week period allowing 
two weeks for a response. A minimum of 25 surveys was needed and it was hoped that 
50 usable surveys could be obtained. A total of 42 surveys were recetved at the end of 
three weeks of which 40 usable surveys were compiled to provide the data set located 
at Appendix B. 

The actual survey used to gather responses is located at Appendix A. The state 
variables, divided into contiguous levels as listed, are organized into the 64 combinations 
or instances. These are then randomized to alleviate any possible influence on the re- 
spondents or judges due to ordering of the variables as each instance is evaluated. The 
number of levels and resulting number of instances has been kept as small as possible 
while still providing sufficient fidelity in the data in order to demonstrate the validity of 
this method without undue effort in the data collection. The cross-country or localized 
speed is broken down into levels of 0 to 5 KMPH, 6 to 15 KMPH, 16 to 30 KMPH, and 
31 to 45 KMPH. The number or degree of LOS is divided into no lines-of-sight. | line- 
of-sight, 2 or 3 lines-of-sight, and 4 or more lines-of-sight. The range is divided into one 
ktlometer bands out to 4000 meters, The exact value used for each variable within each 
band was also determined randomly in the case of speed and range using a pseudo- 
normal random number generator with the median value of the variable in that level as 
the mean. This provides further assurance that judges are not influenced by ordering or 
repeating of values. 

The fixed variables are defined in an operations order given as part of the sur- 
vey. The Central Army Group (CENTAG) area of Germany is the setting for an Army 
‘92 balanced team with M1 Tanks, M2 IFV’s, and attachments to conduct a deliberate 
attack against T80 and T72 Tank, BMP IFV, and BRDM2 ATV equipped Soviet units 
in defensive positions. Information on aircraft, helicopter, and artillery support is also 


provided. 
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The introduction of the survey also contains detailed explanations of the value 
being rated and the rating categories with examples of instances that would normally fall 
fren tnat Category. The five categories are: 


eesO \iovement Potential 


bho 


Marginal Mfovement Potential 


Effective Movement Potential 


ow) 


4. Good Movement Potential 


Excellent Wfovement Potential 


Cn 


This rating of the state of a tactical element with respect to possible tactical 
movement based on the combination of variables or the instance is then combined with 
the corresponding ratings provided by all other judges in order to establish a data set 
from which an interval scale can be determined. This comprises the raw frequency matrix 


located at Appendix B, Section A. 


D. RELATING THE COGNITVE AND PHYSICAL CONTINUUMIS 

The ability to accurately relate the cognitive and physical aspects of the battlefield 
Varies inimensely depending on the action being studied, the scenario, and the target 
population. Numerous methods can be applied allowing for complete translation of the 
cognitive continuum into the physical or structuring the resultant scales in order to 
perform a combinatorial optimization. Two specific methods and the underlving theory 
associated with them are examined in this study. 

1. Computing Interval Scales from Categorical Responses 

The theory involved in deriving interval scales using the Law of Categorical 

duideements 1S fully explained by Professor Glenn Lindsay {Ref. 12: pp. 6-12] and W.S. 
Torgerson [Ref. 10: pp. 207-227]. However, a review of the major points and underlying 
assumptions is worthwhile. There are 6 key assumptions when using this technique. 


1. Judges’ feelings about the scale value of an instance are normallv distributed across 
the population with mean y, and variance o?. 


2. The categories must be a mutuallv exclusive set of successive intervals which col- 
lectively exhaust the property continuum. 


Five categories are required for rating instances. 


4. Judges’ feelings about a category upper bound are normally distributed so that for 
category J, the upper bound would have mean yw, and variance o?. 


s 


This also implies that judges feelings about the distance between an instance value 
and category bound are normally distributed. 


6. Category bounds have the same variances for all } categories. 
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The technique takes the raw frequency matrix for all instances as mentioned and 
first divides these frequencies by the total number of respondents in order to determine 
the relative frequencies for each instance and category. These relative frequencies are 
then summed across the categories for each instance to derive cumulative relative fre- 
quencies which are related to cumulative values under a normal curve. These cumulative 
relative frequencies for category ] of instance 1 become the normal probabilities, p,,. used 
in computing the interval scale. The p,, are grouped based on which instances have cat- 
egories with similar ratings and all p,, greater than 0.98 and less than 0.02 are rolled over 
into the adjacent category in order to reduce the effect and possible influence of outlying 
judges. Treating the p,, as the area under a Normal (0,1) Curve, the Z matrix of normal 
values can be constructed. Row, column. and grand averages are computed from this 
matrix along with an estimate of the standard deviation. The interval scale values of all 
instances are computed as the grand average minus the value of the row average times 
the standard deviation. A linear transform is then used to obtain the correct distance for 
the see: 

The 10 step method of interval scale development given by Lindsay is shown 
below: 


|. Arrange the raw frequency data in tablular form where the rows are instances and 
the columns represent categories. Columns should be sorted in a rank order from 
least to most favorable as left to right columns respectively. 


2. Compute the cumulative relative frequencies for each row and record these values 
in a new table. This table is referred to as the P array and all values of p,, greater 
than 0.98 or less than 0.02 are combined with the following category eliminating 
extreme cases. This creates an n by (m-k) array where k is the number of columns 
removed. 


3. Treating the p,,values as leftward areas under a Normal (0,1) Curve, find the cor- 
responding value of Z from a normial distribution table. Record these values as a 
newetable which 1s the 29 eee) ectuean 


4. For each instance, i, in the Z array; compute the row average, 2— 


For each column or category, j, in the Z array; compute the column average, 0. 
Note that 6, is the value of category j’s upper bound on the interval scale. 


6. Compute the grand average, b, of all Z values. 
m—k 


7. Compute B= ) (b,— 6)? which is the sum of the squared column differences. 
j=1 


m—k 


8. For each row compute 4,= \(z,,—z,)?, the sum of the squared individual differ- 
ences. et 
B ; SIS 
9, For each instance compute “7° an estimate Oi a 


t 
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- a eB . 
iQ) Finally, for each row or instance compute S$, = b — z, x \/ a fore ula, 


ive sya allleseine ie inicrvanscale Valles olathe instances which are on thesame 
scale as the category bounds, 6, . This now enables us to impart distance and a rational 
Origin to the scale through anchoring and linear transformations. The above procedure 
is performed on the raw data at Appendix B. The data is segregated into 7 groups cor- 
responding to how instances were rated by like categories. The first and last groups 
cannot be utilized in determing interval scales since they represent extreme instances in 
the survey, have incomplete Z,; arrays, and only contain one boundry point. All groups 
have incomplete Z,, arrays which is why they are segregated. A tactic which enables us 
to continue with the scale development requires performing the scaling procedure on 
each group and then performing a linear transformation between groups with two or 
more commion boundries [Ref. 12: pp. 18-20]. Each group of data is then normalized 
and an interval scale computed utilizing the A Programming Lanquage (APL) program 
at Appendix B, Section G. The APL program “NORM?” was originally developed by 
Paul Ml. Crawford during his thesis research [Ref. 8: p. 89]. The results of the normal- 
ization and interval scale for each group are located at Appendix B, Section H. 

2. Converting Interval Scales to Usable Forms 

The major requirement at this point is to transform the interval scale values for 
all groups into one contiguous scale by establishing inital boundrv values for one group 
and then performing the linear transformations on the remaining groups. This provides 
Consistent values for all instances evaluated on one usable scale. Another necessity in 
order to use the scale in models is to convert the interval scale into a ratio scale by es- 
tablishing a rational origin and anchor point or points which intersect on the phvsical 
and cognitive scales. This allows us to directly map the two continuums together simul- 
taneously. The boundry conditions apphed will begin with Group II where a rational 
Origin can be applied at the boundry in common with the extreme points shown at Ap- 
memes, section &. A methodsto anchor the other extreme end of the scale will also 


be discussed. 


Ke, 


a. The Natural Origin and Anchoring in Computing Interval Scales 
The extreme cases reflecting excellent potential correspond to cases mem 
reflect little or no anti-potential. We can reverse the rating categories if desired ( which 
was done in this case ) so that excellent TMP or no anti-potential is the worst rating and 
then perform the normalization procedure and develop the scales. This ts done in order 
to correctly place the rational origin. The physical scale of the battlefield based on 
traversal times between cells or nodes reflects only the effects of constraints on the 
physical mobility of a given svstem. At long ranges when there ts no inter-visibility with 
opposing forces and the system can move with great speed, there ts no effect due to the 
cognitive considerations of the battlefield. This zero point or rational origin on the cog- 
nitive scale corresponds to the effect any cognitive rating of these extreme cases would 
have and also the respective value of this boundry on the scale. An arbitrary value of 
100 was selected for the upper boundry of Group II. This now allows for the use of the 
APL program “TRANS” shown at Appendix B, Section G. This program was also ort- 
ginally developed. by Paul Crawierd | REP cas. 70). 
b. The Linear Transformation 
The linear transformation performed by this program solves a set of simul- 


taneous equations on each group of the form: 


Lower Bound = a + ( # x Best Column Average) 


Upper Bound = « + ( # x Worst Column Averag?2) 


where the column averages are already computed and the bounds are input. Once & and 
B are known, the instance values and other intermediate boundry values can be easily 
computed. These other boundry values along with those input are then used to perform 
linear transformations on all other groups with two or more boundries in common. All 
groups can then be combined in order to form one contrguous scale and the value of the 
extreme cases located at the rational origin may be added to the set of scale values of 
the instances. The set of transformed instance values is shown at Appendix B, Section I 
in the last column. The set of instance values along with its corresponding set of vari- 
able values are combined to form the regression matrix shown at Appendix B, Section 
iE 
3. Regression Analvsis of the Interval Scale Values 
Multiple regression techniques are used on the data and derived scale values. 


Two different regresston packages, Grafstat and SAS. are used applying backward elim- 


ination procedures. Plots of the independent variables taken two at a ume with the in- 
stance values of T\[P are done and then surface fitting procedures applied in Grafstat. 
These plots are similar to scatter plots and used in much the same wav. Relationships 
between independent variables are of particular interest. The example plots in Appendix 
C, Section A reflect the linear relationships between pairs of variables. The fitted sur- 
faces appear as planes. Due to this result, multiple linear regression is used to include 
product terms of all pairs of independent variables. Appendix C, Section B gives the re- 
sults of backward elimination using Grafstat and Appendix C, Section C gives the cor- 
responding results of backward elimination using SAS. Each section gives the best 4 and 
5 variable models based on the largest R? with the best acceptable levels of significanee. 

The respective best regression equations for the four variable case and the five 


Variable case are: 
iP = 115824 —0. 01S RNG) + 12.0(L0S) — 2.0.8 PD) + 0.2.LOS x SPD) {1} 


ge = 11232 —olotege\G) 2 11.4(hOS) = 1. (SPD) , 
+ 0.21(LOS x SPD) — 0.00015(RNG x SPD) 12} 


These were the only two equations found acceptable based on the multiple regression 
analvsis. 
4. The Psychometric Power Function 
The “Psychometric Power Function” is the selected best regression equation 
from the regression analysis. The scale can be fully determined by taking the set of 
maximum and minimum values for the independent variables and using these to com- 
pute the corresponding maximum and minimum values of TMP from the regression 
equation. The user must be aware that use of this equation or power function 1s only 
valid for the initial range of values surveyed and used in the regression. The nunimum 
value reflects the rational origin or highest TMP which imparts zero additional time to 
the physical scale as discussed previously. The problem now becomes how to anchor the 
other extreme of the scale which reflects the maximum scale value or no TMP. 
5. Translating Values of the Interval Scale to the Physical Scale 
The essential step proposed bv this methodology requires translating the cogni- 
tive values obtained using psvchometrics into the same scale as the phvsical traversal 
time scale of the battlefield or to be able to somehow perform a combinatorial opti- 
mization between the two scales. There is a fortunate structure in the data resulting 


from the survey. Referring to Appendix B, Section E; there is onlv one extreme instance 


occuring on the end of the scale opposite the rational origin. It is therefore possible to 
prompt the user or survey a selected population to determine the amount of additional 
time in mission completion a combatant would be willing to accept in order to avert such 
a tactical instance. It 1s only possible to construct such an anchor due to the fact that 
there is only one or are a relatively small number of related instances placed at this ex- 
treme end of the continuum. The scale value of this combination of variables or instance 
is easily calculated. It becomes quite simple to then translate the cognitive scale to the 
time scale by dividing the scale or TMP values by the maximum value computed from 
the power function for the extreme instance which results in a (0,1) scale. Then muluplv 
through by the given additional time value from a population mean or as provided by 
an individual respondent. This ume scale obtained using the cognitive data can be used 
addiuvely with the physical traversal time scale. Values for each individual user or a 
mean value for a selected population could be preferred depending on the application. 
This translation allows for speed and simplicity when performing the optimization using 
Dijkstra’s algorithm which is discussed in Chapter [V. Such an approach 1s highly pre- 
ferred if at all possible. Tlowever, it is often not possible to relate the cognitive and 
physical scales as in this research. Other opumuzation techniques which work in combi- 


natonial optimization applications are also discussed in Chapter IV. 


tJ 
tu 


IV. MODEL CONSTRUCTION 


A. GENERAL 

This chapter outlines overall structure of the model and how the methodology al- 
ready discussed is controlled and implemented. There 1s also some detailed examination 
of algorithms and methodologies; particularly concerning line-of-sight, relative effective 
range, and minimum cost optimization. This study follows a logical order beginning with 
required inputs, continues with algorithms and theory used to process data in the model, 
examines optimization techniques, and concludes with the model output. 

Mieeresearch version of this model. called the Tactical Optimum Path or TOP 
Model. is located at TRAC; Monterey, California and at the Mobility Systems Division, 
U.S. Geotechnical Laboratory; Vicksburg, Mississippi. The current model contains ap- 
proximately $0,000 lines of code with most routines done in FORTRAN and the user 
interface, graphics. and data management executed in the C Programming Language 
which is used as a shell. More detailed documentation and a draft users manual for the 
model will be published as a joint technical report through the U.S. Army Waterways 


Experiment Station in late 1989. 


B. TERRAIN DATA BASES 

The digitized terrain data bases used for this model are the identical terrain data 
bases used for the CAMMS (Version 2.0). This is 100 meter grid explicit terrain data 
translated specifically for this version of the CAM MS [Ref. 7]. The model allows for se- 
lection of a geographic region, country. area, and specified “Quad” sheet. This selection 
corresponds to one of the fixed scenario variables, area of operations, identified previ- 
ously. 

A detailed explanation of explicit terrain data and its use can be found in Lecture 
Notes in High Resolution Combat Modelling by James K. Hartman [Ref. 4]. We will 
therefore dispense with any major discussion of these data. It 1s important to note that 
many of the terrain parameters provided in the data set are used in other algorithms. 
This is true of the line-of-sight (LOS) algorithms used. Portions of these algorithms have 


been altered to conform to this specific terrain data base. 


C. REQUIRED INPUTS 

The remaining required inputs primarily relate to the scenario development or fixed 
variables as well. Two of these variables needed as inputs are the appropriate nussions 
of each force. Two other variables comprise the equipment mix available to each force. 
These variables are input to the model by the user during initialization of the program. 
The variables listed up to this point encompass all the fixed scenario variables for the 
methodology. It is important to stress that only the psychometric function determined 
for that specific combination of variables or scenario is valid for that combination in the 
miedel: 

There are presently three other inputs to this model which must be addressed 1n our 
discussion. These are the element locations for each force, current environment, and the 
parameter related to the extreme scale value selection for translating the cognitive scale 
which was discussed in Chapter III, Section D. The necessity of this last input is evident. 
The locations and environment are necessary in determing path start, intermediate, and 
end points as well as range and line-of-sight results. The environment also affects mo- 
bility calculations and may be plaved using real tune data once the Airland Battlefield 
Environment (ALBE) system is incorporated in the CAMMS. 

The user inputs these data to inelude season, current weather, light condition. and 
soil-moisture condition for the environment. The user must input initial locations and 
anv required subsequent locations depending on the mission as each type of combat el- 
ement is entered. The user is then prompted as the last step prior to initiation of the 
program to provide the time value corresponding to the extreme case or anchor cited 
previously or to select the expected value derived from the sample population. This 


concludes all required inputs necessary to run the model. 


D. THE MOBILITY MODEL, CAMMS (VERSION 2.0) 

The CAMMS has evolved from research, testing, and analysis conducted over the 
past 15 years. The CAMMS was initiated in 1983 in order provide a capability for pre- 
dicting mobility in a field environment in real time [Ref. 7]. There are numerous outputs 
available from the CAM MS relating to mobility as listed in Chapter II]. This model will 
utilize two specific output mappings and files from the CAMMS Route Analysis Pre- 
diction: the Vehicle Speed Prediction Map and a version of the Vehicle Time Contour 
Map. These mappings provide the localized cross-country speed or vehicular agility 
which is one of the independent variables necessary for the psychometric function and 


the traversal times along each segment between nodes which comprises the set of arc 


traversal values for the physical continuum. These output values are valid for the speci- 
fied item or Weapon svstem on that terrain under the specified environmental conditions. 
The weather and environment are incorporated in all calculations through the Soil 
Moisture Strength Prediction which is one of the first routines performed by the 
CAMNIS [Ref. 7]. A speed value is stored for all vehicle types for each grid cell or node 
to be used in later calculations. The arc traversal times for segments between all adjacent 
grid cells or nodes are also stored for later use. 

The output provided up to this point by the CAMMS are all mobility modelling 
calculations to be used in the model. The remaining areas we will examine pertain to the 


cognitive considerations affecting movement and performing the optimization. 


E. LINE-OF-SIGHT MODELS 

The model also determines if lines-of-sight exist and to what degree or how many in 
Order tO Obtain the next independent variable, LOS. There are several excellent means 
by which to determine LOS. A detailed explanation of line-of-sight calculations based 
On all tvpes of terrain data can be found in Lecture Notes In High Resolution Combat 
Modelling by James kK. Hartman [Ref. 4]. We will therefore dispense with any major 
discussion of LOS modelling here. 

The LOS algorithms used in this model are derivations of two LOS subroutines, 
IVSCAN and VISTA, developed and used by the BDM Corporation in the Battalion 
Combat Outcome Model (BCOM) and U.S. Army TRAC respectively. The subroutines 
allow for point to point LOS determinations or inter-visibility fan determinations de- 
pending on the desired application and which routine is used. The subroutines are 
modified to use the explicit grid terrain data parameters from the CAMMS terrain data 
base. 

The number of LOS for a given grid cell or node is determined by using the inter- 
visibility fan from each known or suspected threat location to obtain a frequency count 
of the total number of lines-of-sight for each cell, FLOS. This value is also one of the 


independent variables stored for use later in the model. 


F. RELATIVE EFFECTIVE RANGE METHODOLOGIES 

The value of range used in the model requires a more detailed discussion than do the 
other variables. One of the assumptions given initially addressed the fact that the forms 
of the range dependency equations studied and proposed by Seth Bonder and Bob 
Farell are correct |Ref. 13]. Several such forms and their use are expounded bv James 
Taylor [Ref. 14: pp. 93-99]. One similar form was used by Gordon Clark in DYNTACS 


in determing the “Difficulty”, E, from his set of comparative judgements [Ref. 3: p. 98]. 
While there appears to have been a great deal of past interest in such a methodology; 
there has not been anv significant or sufficient analysis of how lethality and particularly 
perceived lethality on the part of a combatant is affected by range. Such a study is well 
bevond the scope of this thesis and 1s not attempted. This is the primary reason the 
stated assumption its made. However, an explanation of the equation forms and why they 
are used in this research is necessary. 

The basic approach concerning range in this model is to assume the combatant has 
sufficient technical knowledge of the systems on the battlefield in order to accuratelv 
assess their potential effectiveness against other systems. It 1s possible to use Probability 
of Hit, Probability of Kull|Hit, and Total Probability of Kill values from classified 
AMSAA JANUS model scenarios along with the classified effective range data for a 
weapon svstem type to fit one of the equation forms previously proposed. This becomes 
the perceived as well as actual lethality curve for that svstem. The systems all have var- 
lous maximum effective ranges with differing curves and equations. The systems are 
equated through the use of a baseline range which is assumed to be the range of the 
battlefield against which the actual range to maximum effective range ratio or difference 
is adjusted. This results in a Relative Effective Range, RER, which reflects the range 
civen in the survey and used in the model. 

Some technical data from unclassified reports, periodicals, and the general popu- 
lation’s knowledge files were easily obtained. These sources provided some idea as to the 
shape and resulting type of equations the lethality curves for various systems possessed. 
Ranges were obtained from current military periodicals [Ref. 15]. The next step is to find 
a set of fictitious vet properly structured and ordered parameters for the equation forms 
identified from these sources. A number of potential forms could be evaluated. Actual 
ranges and parameters exaniuned in JANUS data files or other sources are not used in 


this research. 


The following four equations are examined based on discussions with and previous 


analysis performed by the aforementioned combat modellers: 


RER = RG) eae be 3 
RM, 13) 
RER = RG\e~ x(RM, — eae {J} 
-» _ RG mx RA; _ 
RER= oneal +1 F cos{ Ri {5} 

_, RG Oi | | 
a on +) y+ a RM, -); {6} 


where, 
fees the Kélative Effective Range usedin, the model 
RG is the baseline range used for this tactical scenario 
RA, 1s the actual range to svstem 1 
R.A, is the maximum effective range of svstem 1 


u, a, andy are the respective shape scale parameters for each respective equation form 


Equation number 3 is a basic Power form of the range dependency equation similar to 
the actual versus maximum range ratio attrition equation used in Lanchester models 
(Ref. 14: p.98]. Equation number 4 is an Exponential form which works well when fitting 
non-ballistic weapon systems with fairly constant lethality over range such as current 
generation anti-armor missle systems. Equation 5 is a Half-Cosine form which works 
Well for systems that tend to have much smaller lethalities in comparison to most other 
systems such as smiall arms. Equation 6 is a Full-Cosine form applicable for systems 
which may have minimum as well as maximum effective ranges. [Ref. 14: pp. 96-99] 

Sample graphical results of the plots performed is useful at this point. Figure 1 
shows a comparison of the Soviet IT72 main Battle Tank versus an AT5 Spandrael 
Anti- Tank Missle. 
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Figure 1. Comparison of Relative Effective Range for Soviet NIBT vs. ATGMI 


The tank 1s represented by a Power equation with a uw of 4.0 and RAV of 3200 meters and 
the ATGM is represented by an Exponential equation with an « of 0.00080 and RV of 
5000 meters. The RER is plotted on a logrithmic scale due to the extremely low values 
of RER at shorter ranges. Low RER imphies greater lethality. Note the significant dif- 
ference in Relative Effective Range out to approximately 1200 meters. The tank is per- 
ceived to have much greater lethality at these shorter ranges. No major difference 1s 
evident between 1400 and 2000 meters; but, the lethality of the ATGM remains more 
constant over range and Relative Effective Range changes in favor of the missle from 
2500 meters out to the limit of the battlefield. 

A similar comparison exists for the MIAI Main Battle Tank and the M2,3 
Infantry, Cavalrv Fighting Vehicle where the tank is a ballistic system and the IFV is a 
missle svstem based on the anti-tank missle system in addition to a 25 millimeter chain 


gun. [This comparison is shown in Figure 2. 
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Figure 2. Comparison of Relative Effective Range for U.S. MBT vs. [FV 


The tank is represented by a Power equation with a uw of 4.25 and RAV of 3500 meters 
and the IF V is represented by an Exponential equation with an a of 0.00090 and RV 
of 4000 meters. Note the significant difference in Relative Effective Range signifving the 
tank 1s the greatest perceived threat out to 2000 meters and the IFV equipped with the 
TOW-2 missle 1s obviously greater past approximately 2800 meters. 

Complete listings of the systems included, the equation forms used, and resulting 
parameters are shown in Tables I and 2. Only the Power and Exponential forms are used 
in this analvsis and in the model. Table | lists the values of the maximum range and 
scale parameter, « , for systems represented using the Exponential form equation and 
Table 2 lists the values of the maximum range and shape parameter, uw , for svstems re- 


presented using the Power form equation. 
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It is also important to note that we can compare NATO and Soviet systems using 
these curves. Note that the [72 Nlain Battle Tank has a lower uw value and 1s not as 
effective as the MIA1 Main Battle Tank nor does it possesss as long a maximum effec- 
tive range. The difference ts noticeable when the plots shown in Figures 1 and 2 are 
compared. 

It is also interesting to note how the lethality for various svstems change with re- 
spect to range. This gives some idea as to whether or not the RER calculateduimae. 
sense or is Verifiable. The small change in slope of approximately 0.1 observed within 
the first 1000 meters of actual range for both tank weapons systems is one example. This 
range 1s within the boresight range of each system and there is little or no change in the 
trajectory of tank cannon munitions at such short ranges. The slope and subsequent 
RER support these facts. The slope between 2000 and 3000 meters actual range in- 
creases to approximately 1.7 indicating perceived lethality decreases more quickly. There 
is almost no difference in the RER of 1 at 500 meters range and the RER of I4 at 1000 
meters range for the MIAI Tank. However, there 1s a significant difference in the RER 
of 2600 at 3500 meters for the tank. The TOW-2 missle on the [FV has an RER of only 


1200 at 3500 meters range. While these equations and parameters are fictitious, their 
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application and results correspond to the manner in which experienced military person- 


Mmeltend to compare these systems. 
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For this model, an RER value for each cell or node is determined by taking the ac- 
tual range, RA, for each svstem with line-of-sight to that cell and using this range in the 
appropriate RER equation form for that system. The minimum value of the RER for 
all svstems with line-of-sight, if more than one exists, 1s retained as the RER value for 
eecell. 

The data and equation forms identified and derived from this short study do need 
to be evaluated more thoroughly. However, such an effort is beyond the scope of this 
thesis and would require development of an extensive experiment in order to sufficiently 
determine how perceived lethality is affected by range for various weapons systems. 


Equation forms with exact parameters derived from existing AMSAA data would almost 
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certainly be classified as would the results of such an experiment. This is certainly a 
proposal for future studv and enhancements: but. the assumption stated is used at this 
point in order to continue development of methodology based on the use of 
psychometric techniques without undue requirements due to an analysis of perceived 
lethality and or inclusion of classified data. 

All required inputs to the model are complete once the CAM MS has been run; the 
LOS calculations are performed; the SPD, FLOS, and RER values have been obtained 
for each cell; and the arc traversal times between all adjacent cells are stored. The ap- 
plication of the cognitive portion of the model and transformation of the values to be 


used in the optimization can now be done. 


G. APPLYING THE PSYCHOMETRIC FUNCTION 

The physical time scale is already determined from the CAM MS but we must utilize 
the other data as inputs to the psvchometric or decision function derived from the re- 
gression analysis in order to obtain the cognitive or tactical time scale. The regression 
equation or decision function uses FLOS, RER, and SPD as independent variables in 
order to output a dependent variable, Tactical Movement Potential or TMP. This value 
1s computed as a ratio scale value for each cell through application of a rational origin 
and then translated into the corresponding physical time scale made up of the traversal 
times through the use of an anchor as outlined in Chapter III. Section D. We now have 
two sets of corresponding time values. 

The cognitive or tactical time values represent the additional cost in units of time, 
equivalent to those time units on the physical scale, for a combat element to enter that 
cell based on the existing set or combination of state variables. The two sets of time 
values are additive. Consequently, the tactical time computed for each cell 1s then added 
to the physical traversal times stored for each arc which enters or has a tail node in that 
particular cell. This set of time values becomes the total time along each arc and the set 


of values desired for use in the optimization. 


H. THE OPTIMIZATION 

Optimization of this set of translated network values is performed using Dykstra’s al- 
gorithm which is relatively quick and efficient in this application. However, it is not al- 
Ways possible to determine a rational origin and particularly one or more anchors in 
order to translate the cognitive scale. Dijkstra’s algorithm can solve the optimization 
problem very simply if the requirement is to either reach the objective by minimizing 


tactical difficulty which is equivalent to maximizing PMIP@to reach thes object: ema, 


minimizing only the traversal time. or to reach the objective by minimizing some total 
translated time as in this case. However, the problem is one of combinatorial opti- 
mization when it 1s not possible to determine a relationship between the two scales. We 
will first conduct a quick review of minimum cost optimization using Dikstra’s algo- 
rithm and then examine one method for performing the same optimization combina- 
torially on separate physical and cognitive scales. 

1. Duikstra’s Algorithm 

Dijkstra’s algorithm uses a label setting procedure on a network graph with 
non-negative costs for all arcs in the network. This network structure must be main- 
tained in order to avoid possible negative cycle lengths. Two sets of cells are maintained. 
The first set are those cells whose shortest distance or minimuin cost path is already 
known. The second set consists of all other non-traversed cells. The known set 1s initially 
Semiprsed of only the start or source node. At each iteration. a new cell is added to the 
Known set which has the shortest distance or minimum cost path to the start point. The 
minimum cost value or shortest distance can be recorded or a predecessor array main- 
tained. Eventually, all possible cells on the battlefield are included in the Known set. At 
this point. the minimum cost path from the source to any other cell is available and the 
eeaume can be deternuned. [Ref. 167 pp. 203-208] 

Dijkstra’s algorithm can be made even more efficient bv limiting the nodes 
searched at each iteration. This is applicable in the model since the only permissible arcs 
are those between adjacent cells. The algorithm can also be used to perform a minimum 
cost optimization along a set of intermediate objective points or cells through applica- 
tion of the “Sliding Pattern Concept” [Ref. 5: pp. 22-26]. The optimization is performed 
Ipeallvebetwcen the ordered set of cells n times,svhere n is the number of intermediate 
cells when the order 1s specified and must be maintained. The optimization 1s performed 
in total until all intermediate cells are labeled or in the known set and all are contained 
along one unique path. This can be accomplished using any number of data structures. 
The set of objective points or cells is input by the user when the model is being initial- 
ized. In summation, Diykstra’s algorithm is quite useful for almost any network and can 
be adapted to make it even faster and more efficient with little work. 

2. Combinatorial Optimization 

The optinuzation problem can be solved very simply if the requirement ts to ei- 
ther reach the objective by minimizing some cost along one unique scale when the cog- 
nitive and physical scales can be related as demonstrated in the previous sections. 


However, the problem is one of combinatorial optimization when it is not possible to 
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determine a relationship between the two scales. This requires the analyst to find a 
method by which it is possible to minimize tactical difficulty while also optrmizing the 
traversal time. The optimized time must be less than time T and the tactical difficulty 
must also be minimized within this optimization. Time T should be given by the user 
or obtained from the misston statement of the scenario. The model should select a route 
Which satifies all constraints while solving one equation for both tactical difficulty, rd,. 
and the traversal time, 4. A method that combines both constraints into one equation 


is the use of a Lagrangian multipher, 4 > 0, which can be used in the following equation: 


muin(td + 2 Xx t)x-Ax T {7} 
all x 
wilere. 
x 1S a vector of the arcs which make up a network solution 
td 1s the tactical difficultv based on exposure associated with each x 
fs the travel trme associated with each x 
Ais the Lagrangrain multipher 


T is the maxtmum time limit to traverse the network 


a. Aysumptions 


I. J is larger than the time it takes an element to reach its objective using a shortest 
tume network solution. 


2. The network is undirected. This allows an element at destination B to backtrack 
along its route if that route offers a less difficult route and still meets the /@immmie 
limit to reachithe imal objective: 


b. Nethodology 

The equation is used to calculate the shortest path based on / using 
Dijkstra’s algorithm. The algorithm shown at Appendix I uses a hierarchial adjacency 
list to search for the shortest path between the start point and the intermediate point 
and then from the intermediate point to the final objective. Other data structures could 
be used. 

To calculate the maximum possible value, /,,,,, the largest ¢, from all the 
edges in the network is multiplied by the number of vertices in the network. If /,,, 1s 
used as / in the equation, the objective function minimizes the time traveled throughout 
the network with no regard to tactical difhiculty. When / is zero, the minimum possible 


value of tactical difficulty is found. There exists a 4 such that O</< A,,,, gives the op- 


timum tradeoff in minimizing the tactical difficulty while ensuring the traversal time 
constraint is not violated. 

The best lower bound for the equation is determined by solving for 
maxL(s) which can be solved by doing a binary search from /,,,, tO Ama: AS long as the 
slope of L{/) 1s non-positive then the route selected is a feasible solution with regard to 
T. A binary search is conducted on the interval between /,,, and /,,,, until a specified 
imterval. € . is obtained between /,,,, and /,,,-, [his 4 with e sufficiently small indicates 
the search for the optimum path minimizing tactical difficulty while arriving at the ob- 
jective within 7 traversal time has been accomplished. 

c. The Combinatorial Optimization Algorithm 
A VS-Fortran program was developed to solve the tactical movement 
problem (see Appendix [) and consists of the following steps: 


I. The data set containing the edges and their associated time and difficulty values are 
read into a hierarchial adjacency list. The value of the maximum /, and TDMAA 
in the network are found. 


2. Inititially Dijkstra’s algorithm is called to find the slope, objective function value 
and path for A = 0 and /,,,, = [DMA x (Nodes—1) 


Peeexechieck 1s conducted between /,,,, and /,,,, and if 7,..—4m,n Ee then the program 
is completed and the shortest path with the lowest tactical difliculty has been 
found. 


4. /,,,18 the midpoint between /,,,, and / 


“Max 


sn 


Dijkstra’s algorithm 1s called to compute a slope, objective value and path based 
On ae 


6. If the slope is positive then /,,, becomes /,,,, and Step 3 1s repeated. Otherwise the 
Slope 15 NON-positive, indicating a feasible path and 4,,, becomes the new A,,,,. The 
objective function value of /,,, 1S compared to the previous best objective function 
wotiomoOB ViANe If 4,, obrective valuemis largemthan OBJMAX, /,,, objective 
value becomes the new OBJ MAX and the optimal path, OPTX, is updated. Go to 
prep! 5: 


Once the é value has been reached the optimum path and associated op- 
timal values for time and tactical difficulty are outputted to a file. 


d. Conclusions 


The complexity of Dykstra’s algorithm is O( | V 





*) and the complexity of 
the binary search is O{ log, TD\fAX). The Dijkstra algorithm is utilized with each iter- 
ation of the binary search; therefore, the total complexity of the program is 
OGY log, T DAA) 


To calculate the required number of iterations for a set of nodes, use the 


| . 
Tp t0 determine the interval required and then solve for the number of 
“em ax 


iterations, k, using the equation 


equation - 





_ 5 log, | V | + log, TDMAX+2 logs Tay 
log,2 


{3} 


or use ane of sufficiently small value. 

Changes to the program can be made to decrease the run time if more in- 
formation 1s Known about the nature of the network. In a directed network the 
Dikstra’s algorithm can be changed to stop the search once point B has been found and 
then continue from point B until point C is found. The search could also be restricted 
over ceriain portions of the network based on the network’s structure, greatly reducing 
the number of searches. Another method that could be incorporated is changing the 


binarv search to a slope intercept search over the interval between /,,, and /,,,. These 


min 


methods were not incorporated in this model since it could be used to model several 


tvpes of networks with Varving numbers of intermediate points and structures. 


I. MODEL OUTPUT 

The output for the model is designed to be presented graphically in addition to the 
data files for use in models. Some output is available interactively during model initial- 
iZation. All output is available upon completion of a model run for use in planning, 
analvsis, and to provide an audit trail. The key output is certainly the set of optimum 
path points and total route traversal time. 

Interactive output consists of LOS fans displayed to the screen, a selection of base 
maps which include terrain elevation or vehicle speed maps, element locations, area of 
operation outline, and path point locations. Upon initiation of the model, cells are la- 
belled in accordance with Dijkstra’s algorithm and shown on the screen until completion 
of the routine, at which time the optimum path is displayed. The data are also stored for 


later use. Elimination of the interactive graphics capability when the model is interfaced 
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with a combat model will greatly reduce run time. 

Final output is comprised of a selection of mappings which are the vehicle speed 
map or SPD map, FLOS map. RER map, TMP map, and total time map. The locations. 
path points, area of operations outline, and actual optimum path can be imposed on anv 
one of the selected base maps. Base maps for hard copy output may be selected in scales 
from | to 25,000 up to and including 1 to 250,000. Actual cell values for each state var- 
1able can be accessed by selecting a point interactively with a bit mapping utility. The 


set of optimum path points are easily obtained from a data file. 


ay, 


VY. CONCLUSIONS AND RECOMMENDATIONS 


A. GENERAL 

The working model is resident in the warlab of The TRADOC Analysis Command; 
Monterey, California. [Initial verification runs were conducted on 17 through 25 July 
1989. The model is specifically tailored so as to optimize speed and efliciency on a 
PC-AT 386 computer. A Dell System 325 with VGA-Plus graphics and a high speed, 
wide carriage color printer for hard copy output was purchased to support the model. 
Results and subsequent conclusions and recommendations based on the methodology 


implementation and model runs are addressed in the following two sections. 


B. ~CONELUSIONS 
I. Nlodel Verification 

The ability to perform a validation of the model within the scope of this thesis 
is impractical if validation is possible at all. However. verification in terms of the general 
behavior exhibited by the model and the sensibility of the results must certainly be ex- 
amined. A very strong verification is not even possible without additional support and 
experimentation. Certainly, a more stringent verification or validation should be the 
focus of additional study; but. we must rely on the knowledge of experienced military 
officers and combat modellers in order to accomplish this limited verification. 

There were 16 different variations of fixed state variables initiated for the test 
scenario, each with a different set of beginning path points and objectives. The expected 
value of 5.5 minutes for the scale time value is used in 10 of the model runs. A value of 
,0 minutes is used in 2 of the model runs in order to compare the results to the base 
CAMMS model runs where there is no effect due to TMP. There are no discrepancies 
in these comparisons. There are 4 other mode! runs which use times varving from 2.0 to 
100.0 minutes in order to determine some idea of the sensitivity of the model with respect 
to the scale time value. This is related to allowing the user to input this value based on 
their own risk posture where a large value implies risk aversity. Some output from these 
verification runs are located at Appendix E. A full scale demonstration of the model was 
presented to a selected audience of military officers and professors on 19 July 1989. 
Impressions concerning the model results and and proposed enhancements were solicited 


with an enthusiastic response to the results. 
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2, Output Analysis 

Poem Ode mimic tomuorw> (tremely well in terms of the route selection proc- 
€ss using the expected value and in terms of the sensitivity analysis for varving risk 
postures. The first four pairs of 1 to 25,000 scale plots shown in Appendix E reflect 
movement based strictly on mobilitv versus tactical movement using the expected scale 
femesealue. Ihe last two pairs reflect a change in risk posture where the expected time 
scale value is compared to greater times signifving a more risk averse posture. 

The results are quite evident from these output mappings. The methodology and 
use of a psvchometric decision function does provide a set of path points which appear 
to optimize TMP during movement in addition to mobility considerations. The model 
can also be apphed deterministically, stochastically, or individually using different values 
for the scale time value. These output also appear to support sensitivity of the model to 


Varving risk postures. 


C. RECOMMENDATIONS 
I. Enhancements and Areas for Research 
a. Approved Scenario List 

The scenarios discussed in the development of this methodology are a key 
element in further applications and research concerning the methodology. A complete 
Seno! Scenarios of interest detailing all possible combinations of the fixed state variables 
must be determined in order to structure a data collection effort and target appropriate 
sample populations. There will eventually be a decision function derived for each of the 
identified combinations if the methodology is to be used. Thus effort would be best per- 
ierccantirouen the ciforts of the L.S. Army Center for Armv Lessons Learned (CALL) 
using information from throughout the nulitarv. 

b. Expanded Survey and Regression Capability 

The fidelitv of this survey and corresponding regression Were discussed pre- 
viously. Continuing data collection and analysis efforts should be more detailed in terms 
of categorization and the number of levels or groups used for the independent variables. 
There 1s certainly a balance needed between the amount of fidelity desired in the model 
and the magnitude of the collection and analysis effort. The populations used for these 
efforts should also be expanded and examined closely to insure the appropriate sample 
is surveved. Larger samples would certainly lower the error and increase the power of 
any such test. These actions should be a cooperative effort between the Army Research 
Institute (ARI) and TRADOC Analvsis Command (TRAC). 
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c. Multiple Vehicle Paths and Formation Structure 
The need to model units using a formation structure or by defining a path 
width for each arc on the battlefield which can be used to control movement will be a 
necessary requirement prior to interfacing this model with other combat models. This 
area of research alone 1s practically unbounded. A decision should be made by the Army 
Model Board concerning which approach to use. A formation structure could be applied 
similar to that found in DYNTACS [Ref. 3: pp. 150-155] or developing path parameters 
for the digitized terrain data and then limiting the optimization to only those arcs with 
a munimum acceptable width [Ref. 17: pp. 33-36] are both feasible. It is also possible to 
perform a combinatorial optinuzation using LaGrangian relaxation on three variables 
to include path width. This effort would be best accomplished by the TRADOC Analvsis 
Conunand (TRAC). 
d. Relative Effective Range Determinations 
A sinular tvpe of pyschophysical analvsis needs to be conducted in order to 
determine if combatants’ perceived lethalitv sufficiently corresponds to actual lethality 
curves. The same method of categorical judgements could be used to derive exact re- 
gression equations or lethality curves for desired weapons systems in the model or actual 
classified lethality curves could be instituted if there is strong correlation between the 
perceived and actual results. The data and curves would certainly be classified in either 
case. These actions should also be a cooperative effort between the ARI and TRAC or 
the Army Material Svstems Analvsis Activity (AMSAA) and TRAC. 
e. Symmetric Applications 
There is also a need to determine if psychometric analysis and functions 
done with U.S. and allied sample populations can be extrapolated to and or directly 
applied to Soviet or other forces and if so, at what levels and for which actions. This 
analvsis 1s necessary unless some drastic assumptions are made concerning symmetry 
between opposing force decision processes and tactics in general. This effort will cer- 
tainly be classified and require the assistance of all major intelligence agencies. 
2. Mfethodology Utilization and Future Applications 
a. Tactical Decision Tool 
A major advantage of this methodology used in a model versus a Tactical 
Decision Aid (TDA) such as the CAM MS is the ability of this model to provide or make 
a tactical decision in real time. The CAMMS and other models only provide pertinent 
information to the combat decision maker to assist him in choosing a course of action. 


The model can make tactical decisions in a laboratory or analytical environment and 
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could be used in the field environment given the appropriate terrain data bases and in- 
telligence input. The system and model have been developed to be easily transportable. 
b. Tactical Training Tool 

The model would make an excellent Tactical Training Tool allowing per- 
sonnel to compare their decisions in a laboratory or classroom environment to computer 
simulations based on the model. This application could be performed using a variety of 
mediums. This application could even be done as a field exercise given the flexibility of 
current personal computers. 

c. Combat Model Preprocessor | Subroutine 

The methodology could be used as a pre-processor to provide sets of opti- 
mum path points during the initialization of any number of current combat models and 
then be used interactively as an event driven subroutine to recompute path changes 
during the course of a battle. This is highly desirable in development of svstemuc models 
where we wish to remove the human from the process and decrease variability in the 
battle calculus. The major problems will be to determine which events would drive an 
update in the path evaluation and whether to use the paths at an individual vehicle or 
unit level depending on resolution. 

The model has been designed and structured to make it compatible with 
several existing combat models in terms of the terrain data bases used, the path output, 
and program lanquages used. These include the Battalion Combat Outcome Model 
(BCONI) developed by the BDM Corporation. both versions of JANUS combat models. 
> (OREM, and others. 

d. Research and Analysis Tool 

The methodology has numerous applications and great potential as an 
analvsis tool in current form or as part of a systemic research model. The methodology 
can be applied in evaluating anv type of decision process when the appropriate state 
variables can be determined and a proper sample population identified. Similar tech- 
niques can be used to include these decision functions in models leading to true systemic 
models and expert systems. 

This specific model, once scenarios are determined and analysis performed 
to fully expand the model, can be used to analvze alternative General Defense Plans 
concerning movement. [t may be incorporated into models such as the Obstacle Planner 
System (OPS) to assist in analyzing optimum placement of available engineer resources 
in countermobility or mobilitv roles. These are only a couple of possible research and 


analysis applications. 
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e. The Generalized Value System, GVS 

The methodology also provides a theoretically sound method of determining 
total traversal times for the time domain network of the battlefield which is a necessary 
input to the Generalized Value Svstem. There is additional work to be done in applying 
the methodology to the GVS. This specifically relates to the decision level at which a 
future state decision is made in the GVS versus the level at which paths are computed 
in the model. The problem does not exist at lower levels such as section or platoon. 
However, 1t may not be feasible to extrapolate values of the time domain network at a 
higher level from a number of smaller units comprising a large unit. The use of the path 
width discussed previously may be one solution. However, there would need to be a 
method of determining traversal speeds and other data for aggregated units with com- 


binations of vehicles which the CAMMS and this model can not accomplish. 


D. SUMMARY 

The methodology demonstrates a way to determine optimum, tactical movement 
paths for a specified vehicle and or small unit based on the operator’s cognitive decision 
processes, as well as the physical effects of terrain and environment on mobility. The 
approach uses psychometric techniques inherent to the Generalized Value System (GVS) 
in order to determine a decision function based on the specific tactical scenario and given 
equipment configuration which provides a means to determine the Tactical \fovement 
Potential (TMP) for each terrain cell. This cognitive value in an interval scale can then 
be translated into the same scale as the physical continuum. The cognitive time value 
based on the user’s decision process is then added to the physical traversal times for each 
cell computed from output provided by The Condensed Army Mobility Management 
System (CAMMS). This renders a value mapping which can be optimized using 
Dijkstra’s Algorithm or a combinatorial method could be instituted. The resulting sets 
of path points are optimized for speed time and the cognitive tactical considerations 
evaluated using these psychometric methods thereby rendering a truly optimum, tactical 
movement path and optimum traversal time. The movement path and resulting times 
can be used in a multitude of modelling or analytical applications and particularly in 
determining the time values needed to compute the Situationally Inherent Power (SIP) 
of the GVS. In actuality, this methodology could be applied to almost any tactical de- 
cision process in the development of expert systems and models and shows great poten- 


tial for future use. 


APPENDIX A. A SURVEY OF FACTORS INFLUENCING TACTICAL / 
COMBAT MOVEMENT. 


A. INTRODUCTION. 
1. Purpose. 

The purpose of this survey 1s to obtain an estimate of the degradation to a sys- 
tem’s Tactical Movement Potential (TMP) based on changes in 3 key variables and 
combinations of these variables. The results of analysis to be performed on these survev 
responses will then be used to develop input to a larger tactical decision tool designed 
to provide sets of optimum path points to commanders in the field. Your response to 
each question will reflect the value you place on a given combination of the 3 kev vart- 
Miles With respect to Tactical Viovement Potential (TMP). Your responses as a tactical 
decision maker ’ commander will assist in determining accurate representations of how 
changes in each kev variable affect the rating of its relative importance to tactical 
movement. They may also enable us to deternune how the key variables are correlated. 

2) otructure, 
a. Key Variables. 
There are 3 key variables which have been identified as critical in deternun- 
feeeesvstem Ss 1 VI{P. These 3 key variables are: 


ee sstein Speed based on the cross-country movement rate. 


tw 


. Cover and Concealment based on enemy line-of-sight (LOS). 


3. Range based on the distance to known or suspected enemy positions. 


b. Rating Categories. 
There are 5 rating categories from which to choose in evaluating TMP. 
These are a mutually exclusive set of successive intervals that collectively exhaust all 


possible values of the TMP. These 5 rating categories are: 
1. No Tactical Movement Potential. 
2. Marginal Tactical Movement Potentual. 


Effective Tactical Movement Potential. 


ro) 


4. Good Tactical Movement Potential. 


Excellent Tactical Movement Potential. 


Ca 
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c. Discussion. 

A short discussion to insure complete understanding of the variables and 
the rating category vou should select based on vour tactical experience and judgement 
regarding each scenario is required. We will first examine the key variables and then 
define | discuss the rating categories. 

The 3 kev variables represent the primary factors considered by a vehicle 
Or unit commander when selecting an Axis-of-Advance or Route-of- Movement in oper- 
ations such as an Attack, Movement-to-Contact, Withdrawal, or other operation. In 
each given scenario; which 1s simply a combination of the 3 key variables; the system 
employed will have a varving TMP. For example, a system travelling over rolling, open 
terrain at a range exceeding 4 km. from the nearest opposing force with a formation of 
low lulls preventing enemy observation based on LOS would be capable of moving at 
great speed and in relative safety. This would result in high values for each key variable 
and exhibit “Excellent Tactical Movement Potential. On the other hand waa: 
moving through a peat bog and becoming virtually mired within $00 m. of several Known 
enemy positions with little or no cover and concealment would almost certainly be de- 
Stroved or at least badly damaged and become combat ineffective. This would result in 
low values for each kev variable and suggest the system has “No Tactical \lovement 
Potential”. Varying combinations of the 3 kev variables render scenarios which will be 
rated between these two extremes. This leads to the second portion of our discussion. 

[n order to insure consistent responses between the respondents when 
evaluating each scenario, all respondents must thoroughly understand the definition of 
each rating category. These categories are: 


e “No Tactical Movement Potential” implies the system cannot move through an 
area, a No-Go area such as a swamp. or that the system has practically no proba- 
bility of surviving such a scenario in the judgement of the respondent. This may 
be due to a slow speed. numerous enemy observations at a close range, etc.. Ihis 
represents a tactical situation any prudent commander would never select if possi- 
ble. 


e “Marginal Tactical Movement Potential” implies the system can move at a mar- 
ginal rate and has some probability of surviving the scenario which is low. How- 
ever, this is a situation where although the system can survive; the simple act of 
survival would render the system combat ineffective for a period of time. This re- 
presents a tactical situation which is only acceptable as a last resort for brief peri- 
ods of time. 


e “Effective Tactical Movement Potential” implies the system can move reasonably 
well and | or has a moderate probability of survival and mission success all things 
considered. This represents a tactical situation most commanders would consider 
acceptable in most cases although not preferable. 


e “Good Tactical \fovement Potential” implies the system can move at better than 
average cross-country movement rates and or has a relatively high probability of 
survival and subsequent mission success. This is certainly a preferable tactical sit- 
uation. 


e “Excellent Tactical Movement Potential” implies the system can move at or near 
the maximum cross-country speed with relatively little or no probability of being 
engaged. This 1s the optimum tactical situation. 

3. Instructions. 

The survey consists of three sections. This ‘Introduction’ is the first section fol- 
lowed by a ‘Personal History ° section used to collect population data. The actual tac- 
tical situation and survey compose the ‘Tactical Movement Survey’ section which 1s the 
last but most important section. 

Please insure you feel comfortable with the variables and especially the rating 
categories. Fill in the required personal data as completely as possible and proceed to 
tie survey. 

You will begin by reading a synopsized OPORD giving the Task Organization, 
Situation, and Mission. Once vou fully understand the tactical situation you will con- 
Mime oY alswerneg or rating the TMP for 64 questions © scenarios. Please rate ALL 
scenarios in accordance with vour initial judgement in as timely a manner as possible. 
This 1s done in order to replicate the manner in which such decisions are made in fast- 
paced, tactical operations. Do not change answers once you have recorded your initial 
response unless the change is due to an inadvertent error. Answer ALL scenarios. Thank 


you for your cooperation. 
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B. PERSONAL HISTORY DATA. 
1. Present Rank: 


e lhl C, OO] er dco eee —_ 
© ~ MAS O=4 eeeeee —— 
© CPT, O25 ee - 
® Ll, O-e-0-7 ee a 
6) 23 = EM ua eee tam 
OE Greg iene ee eee — 
Sool OTC OW oe uene es. 


2. Branch or MOS: 
© Armorin@is ain saa 
@ Inlantre see ee eee 
© NOU Cie eee 
eo Alt Deteinse re 
oe ic lC VU e tr eeeeeeres t 
Oo PE NISINGel tere 


@ OUieriws eee 


3. Experience: 


o* Timetommeletive ditt ...c5) see Years Months 
e Time spent in Combat Arms units ..... Nears Months 
e Time attached to Combat Arms units ... Years Months 


4. Opinion: 

Do you feel any other variables should be included in the scenarios which could 
significantly affect the tactical decision logic? 

orcs. 3 NO 


If so, what are these variables and how do they impact tactical movement? 
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TACTICAL MOVEMENT SURVEY. 


1. 


OPORD 1-89, HQ/TF 1-18 XXX 


Task Organization : 


Meant 1amkee 1-15 ANN (-) 
RE eC in.) 2 lane Piet Xt) Team HQ’s 
IpEncinecr coc. wy CY 1 Stinger iat 


SIU tIOn: 


da. 


feet Orces: 

Unidentified forces of the Soviet 111TH MRR are defending tn sector. 
The enemy is estimated to be at 72% strength in men, equipment, and supphies. 
Thev are preparing defensive positions. The enemy 1s equipped with organic 
BNIP IFV’s and BRDM ATV’‘s as well as being reinforced by a TSO Guards 
Tank Battalion and what ts left of the Regimental T72 Tank Battalion. Expect 
Emelmerdiiillem Supper thomune iwsG@ Using 122 SPs. Other artillery support 
is possible. [he enemy mav receive Ilelo support in the form of IIIP or 
EIIND-D units if our operations are successful. The status of enemy fixed-wing 
support 1s not known at this time. The enemy has positions about 5 Kilometers 
away. Some enemy forces have been located and identified {assume vou have 
an Intelligence Overlay). Other larger enemy forees are known to be preparing 
positions behind their security zone. 


Pee mcdis t Orces. 

Team Yankee has been conducting a Movement-to-Contact as part of 
oe and tie JO) pricade. The Scouts and Team X-Ray are to our 
front and have encountered enemy forces in a security zone forward of their 
main defensive belt. Team A-Rav lost 3 vehicles todav. Team Whiskey and 
Team Zulu are on are Ieft and nght respectively. The TF has a Battalion of 155 
See lea a ishe wn Zul Willette ealbreast to secure intermediate 
Objectives Frick and Frack. Team XA-Rav will provide supporting fires from 
their present positions. 


PSSION. 


Team Yankee TF 1-15 XXA will attack NI T 1000 hours passing through 


Team X-Rav vicinity CP66 and between teams Whiskey and Zulu vicinity CP99 in 
order to sieze high ground vicinity Objective Sam. Team Yankee will then defend 
iieplace Until relieved. 


Based on the above given tactical situation; place yourself in the position of a 


vehicle. platoon, or unit commander and give vour rating of the Tactical Movement 


Potential for each of the following scenarios. 
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2. Example Survey Scenarios for Questions 34 through 41. 


34. You are moving cross-country as a part of Team Yankee attacking to seize 
Objective Sam under the following conditions: 


» Range to the nearest known or suspected enemy positron 1s about 1475 meters. 
» Your cross-country speed is less than 5 KMPH or 3 MPH. 
» You are withm the LOS of 1 Known or suspected enemy locatton. 


What do you rate the Tactical Movement Potential of this scenario as - 


( )No TMP (_)Marginal TMP (_)Effective TMP ( )Good TMP (Excellent 


35. You are moving cross-country as a part of Team Yankee attacking to seize 
Objective Sam under the followmg conditions: 


» Range to the nearest known or suspected enemy position 1s 3740 meters 
» Your cross-country speed 1s greater than 45 KMPH or 27 MPH. 
» You are within the LOS of 4 or more known or suspected enemy locations. 


What do vou rate the Tactical Movement Potential of this scenario as - 


(_)No TMP (_)Marginal TMP (_)Effective TMP (_)Good TMP {_ )Excelleniai ims 


36. You are moving cross-country as a part of Team Yankee attacking to seize 
Objective Sam under the following conditions: 


» Range to the nearest known or suspected enemy position 1s about 750 to 770 
meters. 


» Your cross-country speed 1s approximately 20 KMPH or 12.5 MPH. 
» You are not within the LOS of anv known or suspected enemy locations. 


What do vou rate the Tactical Mlovement Potential of this scenario as - 


( )No TMP (_)Maregmnal TMP (_)Eflective TMIP ( )\Good TMP (_ )Esxcellenuiigas 


37, You are moving cross-country as a part of Team Yankee attacking to seize 
Objective Sam under the following conditions: 

=» Range to the nearest known or suspected enemy position 1s 2575 meters 

» Your cross-country speed 1s approximately 42 KMPH or 26 MPH. 

=» You are not within the LOS of anv known or suspected enemy locations. 


What do you rate the Tactical Movement Potential of this scenario as - 


(_)No TMP (_)Marginal TMP (_)Effective TMP ( )Good TP ( )Excellentaiaiie 
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38. You are moving cross-country as a part of Team Yankee attacking to seize 
Objective Sam under the following conditions: 


» Range to the nearest known or suspected enemy position 1s 1440 to 1515 meters. 
= Your cross-country speed is less than 5 KMPHI or 3 MPH. 
» You are within the LOS of 4 or more known or suspected enemy locations. 


What do you rate the Tactical Movement Potential of this scenario as - 


eee) Vareinal PNIP ( jEifective IMP ( j)Good TMP {_jExcellent TMP 


39. You are moving cross-country as a part of Team Yankee attacking to seize 
Objective Sam under the following conditions: 


» Range to the nearest known or suspected enemy position just more than 3500 
mictci 


» Your cross-country speed 1s less than 5 KMPH! or 3 MPH. 
= You are within the LOS of 1 known or suspected enemy location. 


What do you rate the Tactical Movement Potential of this scenario as - 


mow ( )Marcinal (MP ( jEfiective TMP jGood WMP (jExcellent 1 MIP 


40, You are moving cross-country as a part of Team Yankee attacking to seize 
Objective Sam under the following conditions: 


» Range to the nearest Known or suspected enemy position 1s slightly less than 
CoO aMerers. 


weeey OUT Cross-country speed is less than 35544 MPH or 3 IPH. 
» You are within the LOS of 4 or more known or suspected enemy locations. 


What do vou rate the Tactical Movement Potential of this scenario as - 


Om. o PP ( )Mereinal TMP ( jEtlective 7 MP. jGood TMP (_)Excellent TMP 


41. You are moving cross-country as a part of Team Yankee attacking to seize 
Objective Sam under the following conditions: 


» Range to the nearest known or suspected enemy position is between 2280 and 
foe Wenietens: 


me Olle OSS-couminy Speed is approxiniately 400K NIP or 25 MPH. 
» You are within the LOS of 4 or more known or suspected enemy locations. 


What do you rate the Tactical Movement Potential of this scenario as - 


eon ae (oo) viercimal IMP iilective 1 MIR (Good 1 NIP (_)Excellent TMP 
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APPENDIX B. TRANSFORMING THE CATEGORICAL RESPONSES TO 


AN INTERVAL SCALE 
A. RAW FREQUENCIES 
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B. GROUPED FREQUENCIES 
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G. COMPUTER PROGRAMS USED 
THE APL COMPUTER PROGRAM NORM 


VNORMLUIV 

V NORX 

AVNORM3:MTX3NOR: ROWAV;CRAV;:S3;4A3 Bs AAT 3 II 

AVNORM3MTX3;NOR3;ROWAV3:CRAV3:S;:A4A3B3AAL3 LIL 
U<'INPUT THE CUMREL FREQUENCY MATRIX' 

MTX<O 

NOR<NQUAN MTX 

SoM en 

ROWAV<(+/NOR )+(€S[2]) 

ROWAV<(S{[1],1)9ROWAV 

COLAV<(+7/NOR )+(S[1] ) 

GRAV<(4/(+/NOR))= (($[4 )x($021)) 

Se eee VALUES ROW AVERAGE 
oie 
O<NOR , ROWAV 

(j<! 1 


AAT Soo 

II<0O 
i 2 Gal 

AA+(NOR(3: (IIT)1]-,ROWAV )x2 
AAR ee) <A A 
Crist. lo / b2 

ZLAT<(S )pAAT 

AAI<+/AATI 
B<+/((COLAV-GRAV )x2) 
SQR<(B+AAI)*0.5 

Soka esl ii.) oso k 
SSI<GRAV- (ROWAVXxSQR) 
AAT<(S(1],1)pAATI 
GRAV<(ST 11, » 1 )PGRAV 

D<! B 

(j<! pret ge= piie We rae es t 
U<B8 

(j<! t 
U<! fa JE ! 
J<!---------- ! 

D<AAT 

[j<! t 

emcees VALUES = GRAND AVERAGE - (ROW AVERAGE x (BtAlT) & . 
Pe ee ee ee ee ee ee 
Bee ae hE ON See 

out ! 


U<' COLUMN AVERAGES '! 


FPOMONOANE WNP ODOONANHLF WNP OOONA HE WNHOFrPODOON ANE WNP OLILILILILILILILILI 


cme Epes ie Mma 9 Oa ec men come Tl cre I coment pa even 1 ge Jie 1 scm Us We ea Use.) fen Ve Were Je Wom oma ema Lem ile) erent! ole mn een eel en em gem | co cof ane | ann | em Lg ee oe 
OnFr FFE FFE EE FWWWWWWWWWWNHNYNNNNNNNNE PRP RPRPRPRPRPRRPOON OW FEWNEH 
PgR to CSS Ses eet La LI Lana Serre eS eee rere ic rcs 


VNQUANCOIV 
V Z<NQUAN P:A:B:C:D 

C1] >((+/(A<(P<0)v (P21)))>0)/L1 

[2] C+ 2.515517 0.802853 0.010328 

[3] D+ 1.432788 0.189269 0.001308 

Cu] P<+((A<(PS0.5))xP)+((P>0.5)x(1-P)) 

[5] B<(@Px 2)*0.5 

[6] Z<+((2xA)-1)x-B-((Bo.x 0 1 2)4+.xC)eC1+((Bo.x 1 2 3)+.xDdD)) 
+0 


Meee b1:<'THERE IS NO QUANTILE FOR P = ',8A/P 
V 


THE APL COMPUTER PROGRAM TRANS 


YVPAZNSCOUIV 
VETRANS 
i) AVTRANS: UP;LOW;:BETA: ALPHA3;MX;COLUP;COLLOW:TRCOL 
eZ H<+'INPUT COL AVERAGES' 
[3] COL<O 
C4] S=pCOL 
Pe ) COLUP<COL([S] 
C6] COLLOW<COL{(1]) 
(7) O<'ZINPUT THE VECTOR OF BOUNDRY VALUES TO BE TRANSFORMED ' 
[8] MX<Q 
ES | O<+'INPUT UPPER BOUNDRY VALUE! 


E10 ) UP<O 
ii | O<+'INPUT LOWER BOUNDRY VALUE! 
[12] LoVv<O0 
C13] BETA<(UP-LOW )+ (COLUP-COLLOW ) 
Ci4) ALPHA<UP-BETAXxCOLUP 
fas } TR<ALPHAt+EBETAXMX 
C16] TRCOL<ALPHA+BETAXCOL 
is 7 | O<'TRANSFORMED UPPER BOUNDS ARE'! 
[18] C<Z7RCOL 
[19] O<' ! 
p20 } U<+! TRANSFORMED DATA ! 
] 


U<TR 
V 
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H. THE NORMALIZED AND TRANSFORMED VALUES OF THE INSTANCES 


TCROUP I 
NORMALIZED VALUES 


~1.43980047 


~0.7551784916 


~1.150435626 


~1,.64521144 
71.64520004 
~1.43980047 
o1 . 64S 
51. 6U Sas 


~1,036431485 


ROW AVERAGE 


_0.252Gne2672 
~0.5240018704 


0.5240018704 
0O.8414567174 
Opes 1eceo2s 
Ove tet essio7s 


SO -125sc0eSo7 
70.384877085 

70. 887550204 
Molise 7/ sis oto 
~O.597404U8985 
~0.06254483635 
"0.2529332678 


COLUMN AVERAGES 


“2.198129 767 ©2101829137 727-1 17 24 es 


GRAND AVERAGE 


P.43 5567 seo 
2.40935 2es25 
ee 2 Lo olor 
2,891607 505 
2.17 36ei oc 
2.583007 09 e 
3, OOU Ste a2s 
2.6101949UuU8 
2.098 1 50acs 
S205 noe LOG 
1.779340999 


SCALE VALUES 


GRAND AVERAGE - 


1.43980047 
1.64521144 
1.43980047 
1.64521144 
O.. Satara 
Orr so 1764S ike 
1.036431485 
O.8414567174 
0.67418914 
1. OSG Skt Ss 5 
0O.8414567174 


0.5702896907 
0.6542 2208a7 
0.1060666208 
0O.40274H2703 


0.11377 tee 
~O0,.42uU970G Fr 
~0.2658453318 
~0.2623665976 
~0.5228090662 
"0.2297 7482s 
~0.1493026786 


(ROW AVERAGE x (B+AI) x .5) 


mO-c00se5u7 53 
~0.7117889701 
~0.1054995408 


~O.4087904H96 


Oi DOS a eo ay 
0.389057H987 
O. 2082 S417 76 
O24 a2 eos 
0. 5004515703 
OIA ARS repens ee) 226! 
O2.16 2456339793 


~0.02086518298 
0.02086 5a62c8 
~0.02086518298 
~0.02086518298 
~0.02086518298 
~0.02086518298 
(0, 020865128208 
0.020e6S 18295 
~0.02086518298 
~0.02086518298 
~0.02086518298 


Ofe7 OZes boy 
Oj 650227 0257 
0O.1060666208 
O40 2782 703 


20. ts 7 ke 7 Yt 7 
HO. Woue 70012 
~0.2658453318 
~0.,2623665976 
~0.5228090662 
SO. 220 7 Aue sou 
~0.1493026786 


1.366884769 
1.056099744 
0.7979356aiGe 
0.9632049299 
1.110947 s5e2 
0.9645920205 
0.862/1727Ge 
1.013807ee 
0.9971453S37a 
0.85991400964 
1.22/8986e2 


COLUMN AVERAGES 


Memos 129767 0.01829137737 1.1172428H1 


PRANS 


INPUT UPPER BOUNDRY VALUE 


100 


ae 
INPUT LOWER BOUNDRY VALUE 


0 


TRANSFORMED UPPER BOUNDS ARE 


O70 


meeeeot27358E 15 53.53673384 100 


LPRANSFORMED DATA 


17.17841399 
21.00486097 
feet O2545 

eee O91 24367 
Se, s04U7038 
eee 6 997163 
ee 75108343 
geo SIONS 

fee 25103621 
mere 14895057 
98.76346389 


Benogrl 


NORMALIZED VALUES 


0 
~0 
0 
me 
ml 
78 
~0 


1887560404 
~384877085 

mweo29332678 
melo Lags625 
mr oO0E35626 

72240018704 
mo LSIl998625 


1.43980047 

ico l/ 26707 
i koO4S S626 
i OSG silts. 


0.9345033954 
O./55 17 SSG 
O.7551784916 


COLUHN AVERAGES 


~“O.4uU82005164 1.0504652uU5 


GRAND AVERAGE 


fee 26098153 
me OO 8 O15 


AVERAGE 


0.2 Soe 2G 
~4YU8U258358 
24487511792 
56) ons) diese ales 
see? OGiGd 1.5 3 
eos ssesi06 
~2184893146 


O96 7 22002765 
0.917513044 
Valeo i is) Selle, 
0. 6dS4d511 332 
O”,. 97607 0sts a 


SCALE VALUES = GRAND AVERAGE - (ROW AVERAGE x (B+AI) x .5) 


"O42 73oS 25 0.3011323644 0.6255222146 
~0.1021066199 0.3011323644 0O.4484258358 
~0.1780916163 0.3011323644 O” tl Soa ae 
~0.09616729347 0.3011323644 0.3591158114 

0.3787390085 0.3011323644 ~0.1079661155 

0.1657110832 0.3011323644 0.1155883106 
~0.003925450303 8.301 tezceue 0.2184893146 


COLUMN AVERAGES 
~0.4482005164 1.0504652H5 
TRANS 
INPUT UPPER BOUNDRY VALUE 


Ge 

BS eooe7 24 
INPUT LOWER BOUNDRY VALUE 
Els 


Or0 
PRANSFORMED UPPER BOUNDS ARE 
Se o527 e798 15 53. 536764 


LRANSFORMED DATA 
Bazvolte SZ 

12 .26eoo48 
5.645908 1i07 5 
2072 bees 
PSPs ON 10) Se. 
2l.9307 2188 
5,07 OSGeo 


EGROGL LIT 
NORMALIZED VALUES ROW AVERAGE 

Bier pe dc abe ssi 2/5 1.150 356.26 O. Ut lodsigee19 
~0.67418914 O.8414567174 0.082624 7S5866 
~1.036431485 O0.5974048985 —0.2195132933 
~1,036431485 0.67418914 "0. dee a7 2G 

0.1887560404 13281722737 OU oer eS 5 
~0.5974048985 O. 7532765 ste 0.07888679654 
~0.67418914 0.1887560404 ~—0.2427165498 
~0.3181998625 1.64521144 OO”, 6605057688 
(0.452333367.0 1.1505356265 OF dt ss oui 23 Z 


~0.5974048985 OO. 7507 ee sas 0.07888679654 
~0.8414567174 0.4533333679 0.1940616747 
7037 921763516 0.5240018704 ~0.1155883106 
~0.9s4 5033954 1.43980047 0. 2Zo26465 0 / 1 
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0.9202417922 
0.8992322745 
1.067905 7 26 
1.10632 73233 
0.7183 03se0> 
1,.171582384 

1.396203 2e3a3 


71,.036431485 O. cuismeasise25 i076 590158114 


~1.64521144 VO, eUetssee25 Vee 0i56 513 
_1.43980047 0.5240018704 ~0.4578992996 
1.64521144 Oe Uceottecec5 TOry cecsg019 


COLUMN AVERAGES 


Momei 26365611 0.7084504uuN96 


GRAND AVERAGE 


7,.078445099 
1.148591183 
ijeeost710664 
1.463111462 
0.597294679S1 
0.9147409136 
Ors 723371922 
1.927491972 
1.286037493 
Srrord74uO91s6 
0.8382406825 
oo 181511992 
2.818659422 
217513044 
0.8804798637 
Weeo2s259815 
fees 821575 


SCALE VALUES = GRAND AVERAGE - (ROW AVERAGE x (B+AI) x .5) 
Mee 650723627 §6©0.05209305576 0.4161178819 1.0385714u459 
Memrcs60270889 0.05209305576 0.08363378866 1.00358999 

Beeeo 22718504 §0.05209905576 ~0.2155132933 Oreo 10257 
mee sS6022%5 i 0,05209305576 0.1811211726 0.8892018419 
ae 75327804 BOOS 20980 eo76, 027352502085 dee selene) 7 
~0.1408074984 —0.05209305576 0.07888679654 1.124579099 
Oee757359646G 0.05209305576 0.2427165498 iewiezooos 39 
Bemeso61218971  0.05209305576 0.6635057888 0.774716438 
Memos26747026 0.05209305576 0.3485511292 O.9484452039 
~0.1408074984 ~0.05209305576 0.07888679654 1.124579099 
Me 56857464 0.052093055765 0.19406167H7 1.1747 74a 
ees >eese7i12 0.05209805576 0.1155883106 Ie alesS) 1 Jl (ils, 
Mere 109511961  0.05209305576 0.25264385371 O.64064U5468 
mmeeri505256 0.05209305576 0.3591158114 ime 287 8559 
1.07318704 BOeos 205205576 POeceiyese5i3 1.-186249992 
0.3025784214 §0.05209305576 0.4578992996 Oa77t5e21745 
Meeeree 7421652 §~0.05209305576 0.79133393019 Gee sWaIs7622 
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COLUMN AVERAGES 


~0.8126365611 0.7084504u96 


TRANS 


INPUT UPPER BOUNDRY VALUE 
Be 


LOO me 


INPUT LOWER BOUNDRY VALUE 
(ie 


aol e Sileh/eo it! 


TRANSFORMED UPPER BOUNDS ARE 
935 3360 3 ene 


TRANSFORMED DATA 


G35 
To. 
oye! 
7GIS 7S ooo 
us. 
74, 
ooh. 
0 Cle ca oer 
66. 
74, 


81 


Sil 


85 


8Q. 
(ele 
06 Seca 


89 


ie. 
Bix 
cliche 


SCsG3262 
20451078 
OAC ere 2 


51254274 
05848734 
§3687224u 


67 Oars sez 
05848734 
(322266 
96684485 
S2t Za 2 


141278 
SO Ze aUs 
IRS 3 5) 0 19) Ae 


BET<(100-53.536734 )#(1.3298354u4u1-.06802605834u) 


jena sl 


36. 82072982 
ALP<100+(BETx.06802605834 ) 


Ae 


102.5049052 
UPE<ALP+(BETx1.267100775) 


UPE 


1USetesOte, 


TGROUPIV 
NORMALIZED VALUES 


ROW AVERAGE 


TO 8terus6 7174 
"Oe sey s4Si6 


~1,4u3980047 
Ty Stee Ss Be OL 7 


0.2529332678 1.64521144 

0.384877085 1 usseuou, 

0.1887560404 1.43980047 
~0.2529332678 1.150435626 
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72. @) 2 


2.302229 ape 
3 5S64U9DSG877 
-06291868014 
~l80 7660s 72 


Belo 591144 


Mele 281728757 
et O0435626 


mao 452114 u 
me, 63980047 
mit 3980047 
mee 1 39800uU7 
mee 4 39800U7 


~0.3181998625 


Os 252913226 7.8 


~0.3181998625 
ON Seal eras 
~0.06254u483635 
~0.3181998625 


OVS tS Tee sie 25 


~0.1887560404 


COLUMN AVERAGES 


Bo. OSsO02605c3471.267100775 


Bie 3 29835u4uu1 


GRAND AVERAGE 


eee OG 5248825 
2.410174235 
#- 169803347 
meeOo24 7368 

3.911164149 
eee 23920871 
ieee os02612 
Zee? 1LOUS262 
4.7768543 54 
pez 23551553 
u.789607793 
pmo 29996505 


SCALE VALUES 


GRAND 


AVERAGE - 


Paefaor3 56:26 
1.43980047 
0.67418914 
0.384877085 
1.64521144 
1.43980047 
1.64521144 
dr oO sa O26 


BO.7709918921 


WF 13./ VOL6G0 2 


~0.2648154u96 
~0.6718376156 


0O.04762204472 


~0Q0.1060666208 


O.17453694U3 


 Omlisee425279 


(ROW AVERAGE x (B+AI) x .5) 


0O.4106074699 
eer onSsZ211892 
e. 100175096 

We 13/4601038 
we 2080690671 
~0.1739794959 


Mea32sst3213 


Oreos /927 532 
~0.08360358382 
freootg1141717 
~0.190054uu3a98 
071161893582 


COLUMN AVERAGES 


mie 29825041 
LRANS 


~0.04358690797 
BOWOUSSSo907 ey 
ono 4 35 S690 747 
~0.04358690797 
Somos Soo ser 27 
~0.04358690797 
~0.04358690797 
~0.04358690797 
Bo ,043586907¢7 
~0.04358690797 
~0.04358690797 
nO. 043536510 797 


Ooo 22255302 
0.3564996877 


0.06291868014 
Ont 207550371 
~0.2709918921 


O-13700766G¢ 


~0.2648154u96 
~0.6718376156 
O2OUF5 2204472 
~0.1060666208 


O.1745369443 


On 1 593736279 


Bom SO 26058241267 1007 75 


1.041993186 
ielez7so 208 
SO. S993861263 
tO Ode sia 20 
OeI256'37617 2 
0 .. Some 591802 
IeiZzileseZ 

Tee /el17 OSes 
O238402973052 
0.8947048978 
Oressl/7sies 
OO Zo 7 2 


INPUT UPPER BOUNDRY VALUE 


Be 

UPE 
INPUT LOWER BOUNDRY VALUE 
Be 


Dow so? o4 
TRANSFORMED UPPER BOUNDS ARE 
53.936734 99.99999999 749.1630147 


PRANSFORMED DATA 
el fees Cita a) 2 I 2S 
Src Ose 72 
98.81618467 

107 -3665614 

LTO mooe 7S 2 
96.09850519 

114.7626146 

132.4709759 
99.42639298 

104 .3943316 
Co oie dic 7 

TOG. 7833405 


TGROUPV 
NORMALIZED VALUES ROW AVERAGE 

1.010066757E 7 1. 6H S240 20 8.226057706H 1 
56.7426 ore008 > 1 1.439800470E0 3.82805664U8EF 1 
“4. 253 000Sa78 71 1.645211440F0 7.99915223 520 
~2,.5293392678e8 1 1.645211440F0 6.961390862EF 1 
~Q9,345033954UF 1 1.439800470£E0 2.526653 71d 
75.9790 Meee Se 11 765517 SU STBer a1 7.888679654E 2 
~1.036431u4U85EF0 5. 2400S 7 OUR 1 9 2256720 S07 eee 
SC, Uiise717UE 1 1.150435626E0 1.54UUSSUSUSE 1 
SS U1 sey 7 4H ed 1 65211 OO 4U.018773614E 1 
=1 , 036 gates o 2.529332678E 1 ~3.917491087E 1 
~4.150435626E0 1.150435626E0 0.000000000E0 
~4.150435626E0 1.8875604U0UE 1 Petee0e3s97929F"1 
7 121 2S 75 7 EO 7.551784U916E 1 ~2.632751325E 1 
E1.64527 17454080 2.529332678E 1 ~6.961390862E 1 


COLUMN AVERAGES 


_—__ enw es wes Se we we ee we we we we wena = = 


"On 82622656749 1.035027 324 


GRAND AVERAGE 


12732270255 


70 


meso 3560175 
@e2stt/6035 
io 7498782 
1.801476666 
mes LO659422 
Wei 7409136 
me217476129 
Meso sc1/755u4 
om 1759263 

ees 12307332 
2.64700426 

treo oo/171601 
2.074495569 

1.801476666 


SCALE VALUES 


Mes 26285789 
mo 232681524 


al 
7 


~O.6044u809054 


Moo 73261130 


~0.09369293795 
~0.004189260708 


2 


0O.4099845144 


70.039992517 
0.196443159 


S 
st 


0.6698908224 
0.104%3678548 
fe 7 26743588 
O.,.344 9462429 
wee 3 9968397 


GRAND AVERAGE - 


O. 104 ea 7e5ts 
O.1043678548 
0.1043678548 
0.1043678548 
OOS soyco4s 
0.10 36772548 
O, tO4SE 7eots 
OOH S6y-ea 
O.1043678548 
O. 1046672548 
0.1043678548 
Oj. 1O4S67 ess 
O.1O0U 267 S548 
O. 1O0N Sey S54 6 


Ovee 26057706 
O2eeze0S664s 
Ones 2235 
OS esipila see 2 
Omer eteos7 1 


O20 7sose/7/9654 
~0.2562148074 


O.1S44S945US 
Gyre 7.7 36.14 


~0.3917491087 
0 

70.4808397929 
G. D6 A/a AG 

~0.6961390862 


COLUMN AVERAGES 


Mo. e262856749 1.0350213984 


TRANS 


INPUT UPPER BOUNDRY VALUE 


O 
OO 30 


Cre 
INPUT LOWER BOUNDRY VALUE 


LPRANSFORMED UPPER BOUNDS ARE 


meee 14S .1630147 


TRANSFORMED DATA 


99.99999698 


iio /89657 
P@e.5813814 
ier. 2911138 
ios 200945 
pees 141701 
el 38113 
i20,7684663 


7\ 


(ROW AVERAGE x (BtAI) * .5) 


ieee oi sees 0S 
On ecot7 1243 

ioctl Ze Sg 

WR e0Ss2 e135 
OR eodse0152 
diye GIS) edb) eye 

mis 23i4isy 
0O.9344415953 
O.74¥85144S4Y 
1.443584567 
OFoWsgo/7 5013 
Pacey solc 
Gro s/h? 7.6.2 
Ol ithe 8 Seb sya cpe 


116. 6361497 
pees) 2) Alis: 7.) 0 il 
1.24) 736405 Or 3 
142 233343 
130.9959445 
142,60 i250 03 


~J 


tu 


]. THE REGRESSION MATRIX 


QUES RANGE LOS x SPEED xX IMP 
NUN (METERS ) MOCO KOT Y, 100 (KMPH ) RATING 
ol 58) 5/0 0 2600 an 20 See) 
i 3400 1000 2400 P2yoostecte 
ils: 3450 1G ONS) 6) 73//0)/8 9.64908167 
21 2600 0 2500 go 5) 7) Bap ers18) 6! 
54 ILGHS a) 0 2400 293.54070381 
As! 25.5.0 lee 4UOO0O0 Ole S507. 21 Sic 
26 1480 0 3900 1415) (18 13 Oens) 
20 Ove 0 4200 I len hey eps: 
3 3200 2500 4UOO0O0 Zi Oou seus 
36 yoo 0 2000 47.1902545 
as gigs © 0 11 0:0 3409124367 
UU See 1000 4200 Be. out Ose 
ay g50 0 Zo00 1200 Gis oS 37 es 
U6 eas 1000 4000 S(O Fee) dO s Sihhs 
U8 500 0 1000 SZ. sco 2ued 
Jal 3450 250.0 ZO © TSO O 7G 24 
35 3740 4U500 U500 59.14895057 
Zo 3600 0 500 BSo.7estosss 
10 1500 1000 2500 Go wO 305251 
U3 3.300 1000 1000 74,20451018 
18 1450 2500 2600 Som Ot OS 7632 
27 5200 4500 208 Saoe sy Sioa 5 
Ug 2480 0 A ON 4U5.51254274 
a7 Zo 7xay() (0) 4200 74,0S84U87354 
39 3) 3).o18 1000 5:00 89.8368722u 
42 Jb 2 © 0 1200 Gen ibo7 Sao 7 
a6 2600 1000 2200 6On070s 69512 
al 2400 Zo00 2200 74.05848734 
293 2600 0 SIG Sap S22 288 
2 1600 0 500 80.936684UU85 
ue? 2450 1000 1000 PA ae 2s 2 Age 
ud 2305 4500 4000 SS.08 583313 
0 dh 00 4500 4000 a er Zee 
38 rae, 1000 al 2010, S7,oUZiIS04S 
SiG: 2640 2500) 1100 He .23ess0 Ls 
24 so 2906 4U200 Barooe 2a) 20 
64 Spe 1000 2200 Ss 0 O08 72 
U7 3500 4500 1200 Ss.clols4o7 
20 st? 5 2500 20/16, A107 (5) Sy shel ses GL: 
iS 3400 4500 S00 ADAG 6 Wee nelg) (S) 
3 2500 LOIS 500 96.09850519 
56 500 4500 4000 Pie oe 26106 
a) 400 4500 23.010 care eh Wok” 
Zo ZZ 0 USO00 2200 99.42639298 
8 640 2500 2400 HOU. 3 S433 16 
eZ 800 0 S00 Som o 065618 7 
28 660 1000 1000 Oba 7633145 
BZ 1500 25010 2400 99.99999698 


1 


34 1475 1000 
22 £525 2508 
14 TS ss 2500 
60 540 1000 
3) 7 700 2200 
U5 232 2500 
30 LSS 4u500 
16 700 4500 
6 1400 4500 
L P20 2500 
61 PLS V0, 4500 
Shes) 1480 4500 


115.676 Soy 
10's .S8icieas 
106.5517 iG 
1719 .350033) 
121.714 7768 
132.6526 0 
120.7684663 
116.6362 
139.51876am 
124. S825 ere 
142. 23oa Se 
130.93594uU5 
142.612o002 


APPENDIX C. REGRESSION ANALYSIS OF COGNITIVE FACTORS 
AFFECTING MOVEMENT 


A. SURFACE PLOTS 

Surface plots of all possible combinations of independent variables were done in 
order to determine if there were any particular relationships between variables with re- 
spect to TMP which needed to be accounted for in the regression analvsis. These are 
used in a fashion similar to scatter plots in simple regression. Note the linear or planar 


relationship between most combinations shown. 


5—D PLOT OF IMP VEKSUS KANGEANDRS oom 
CONSTANT LINE OF SIGHT = 0 
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Figure 3. Surface Plot of TMP Versus SPD and RNG for No LOS 
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5—-D PLOT OF TMP VERSUS RANGE AND SPEED 
CONSTANT LINE OF SIGHT = 2 OR 3 
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TMP IN INTERVAL SCALE 


Figure 4. Surface Plot of TMP Versus SPD and RNG for 2-3 LOS 
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3—D PLOTOR MMP VERSUS KANCESANS SPEED 
CONSTANT LINE OF SIGHT = 4 OR MORE 
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TMP IN INTERVAL SCALE 


Figure 5. Surface Plot of TMP Versus SPD and RNG for 4- LOS 


3—-D PLOT OF TMP VERSUS RANGE AND SPEED 
ALL LINES OF SIGHT 


















CRUBBRRTIN 
|g pL 
SOE SOT ADOT et 
RETA SELBRRNIN 

LE NSESSAE SSS OS) a 












TMP IN INTERVAL SCALE 
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Figure 6. Surface Plot of TMP Versus SPD and RNG for All LOS 
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IMP VERSUS SPEED ANG NUMEEX OR BINes OF Signi 
ALL RANGES 


TMP JN INTERVAL SCALE 





Figure 7. Surface Plot of TMP Versus SPD and LOS for All RNG 


TMP VERSUS RANGE AND NUMBER OF LINES OF SIGHT 
ALL RANGES 


TMP IN INTERVAL SCALE 





Figure 8. Surface Plot of TMP Versus RNG and LOS for All SPD 
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B. REGRESSION USING GRAFSTAT 
MULTIPLE REGRESSION USING BACKWARD ELIMINATION 
TABLE CRVCCEFPPIC IE ie 
63 OBSERVATIONS &R-SQUARED = 0.9048 STANDARD ERROR = 12.486 


5 VARIABLES ADJ R-SQUARED = 0.8983 

COEF ESTIMATE STD ERR RESTA? SIG LEVEE 
INTERCEPT Peusgu es? 5.2563E0 27 .uU95 1.6653E-16 
RNG3 -1 U7UUE-2 1.4786E-3 -9,.9715 3.4722E-14 
LOS3 1.1545E-2 1.6888E-3 6.8366 5.5107E-9 
SPD3 -1.9871E-2 1.7745E-3 -11.198 6.3838E-16 
(LOS3-SPD3) 1.9839EF-6 6.8170E-7 2 alae 5.1154E-3 


The assumpuons of normality in the residuals and homogeneity of variance in the 
residuals are supported. The model 1s valid based on the standard error of the estimate 


and the extremely highs 


ANOVA RESULTING FROM THE REGRESSION 
ANALYSIS OF NORMAL DISTRIBUTION FIT 
DATA >: RESIDUAL 


SELECTION >: ALL 
X AXIS LABEL: RESIDUAL 


SAMPLE SIZE : 63 
CENSORING > NONE 
FREQUENCIES : 1 


EST. METHOD +: MAXIMUM LIKELIHOOD 
CONF METHOD : ASYMPTOTIC NORMAL APPROXIMATION 


CONF. INTERVALS COVARIANCE MATRIX OF 


(95 PERCENT ) PARAMETER ESTIMATES 
PARAMETER ESTIMATE | LOWER URBET MU SIGMA 
MU 2. OS MO AG 2.9989 PL Bs| cite) 2 2.2781 oe 
SIGHA a 2 Sie 1 “Pe - 5-0 a -eeee O  O a2 0 dy ee 
SAMPLE _ Fai DIG 1 Gs Sy OY GOODNESS OF FIT 
MEAN : ee a ee edn 20 Oe eee 
rs VIR 8 chil! a 1.207 6ee ds oso CHI-SQUARE : 2.9566 
SKEWNESS: -3. 16645 -2 0.0000E0 DEG FREED: | 
KURTOSIS: 2.512020 a0 UO ag SIGNIF >: O.,3oeae 
KOLM-SMIRN : 0.060943 
PERCENTILES SAMPLE BI TTED SLCvi kh >; 0.97341 
ay Se Or OD =I. 37 0054 CRAMER-V M : 0. Gis 7ezieee 
10: “14.305 poe es) So) ial SIGNIF : > de 
Pa ah =8.0305 =6 07 GeL0 ANDER-DARL :; 0.22408 
20 = iO Oe dle Zea ee SIGNIF : > wie 
(se 10522 8.0768E0 
90: 15.320 ble Ceol sia) Jt 
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Figure 9. Residual Analysis of the Four Variable Regression 
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GOODNESS OF FIT ON THE RESIDUALS FROM THE REGRESSION 


CHI-SQUARE GOODNESS OF FIT TABLE 


LOWER UPPER OBS EXP _ eee (CO=E )x20ee 
eg ce iS. 632 6 6.0403 0.040285 0.00026868 
io. S 39 7.8194 11 10.149 Ons 202 0.071344 
7.8194 0 16 deers oil at _ Ome e a7 0.037023 
0 eo gt, 11 tbsp Clb AL 4.3107 TZ Anai6 
7,8194 Ws 16,8) SI 14 10.149 3.38508 1.4612 
15.639 +INF. S 6.0403 1.0403 0.179a6 
TOTAL ae G3 2.9560 
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Figure 10. Residual Analysis of the Five Variable Regression 


C. REGRESSION USING SAS 


1. THE FIRST SAS MULTIPLE REGRESSION 


OPTIONS LINESIZE=80;3 
DATA ONE; 
CMS FILEDEF DD1i DISK GDATA22 DATA A; 
INFILE DD1; 
INPUT RANG LOSS SPED TAMP; 
PROC STEPWISE DATA=ONE;$ 
MODEL TAMP = RANG LOSS SPED / BACKWARD MAXE&; 





SAS 6-50 TUESDAY. PONE. 13, 1388 
BACKWARD ELIMINATION PROCEDURE FOR DEPENDENT VARIABLE TAM. 
STEP 0 ALL VARIABLES ENTERED R SQUARE = 0.89099524 , 
CEP) c= u, OOO OCGine® 
DF SUM OF SQUARES MEAN SQUARE F PROB>F 
| 
REGRESSION 3 84U699.83472023 28233.27824008 160.75 0.0001 
ERROR 59 10362.21574916 175.63077541 
TOTAL 62 95062.05046940 
B VALUE STD ERROR PRYPE I) ss F PROB>F 
INTCPE Gio? 35145 
RANG Seo eel og 0.00156910 15677.32260814 89.26 0.000% 
LOSS 0.01563357 0,.00099485 43371.36043090 246.95 0.0001 
SPED -0.01595209 0.00122666 29702.31143972 tc uam 0.0004 
BOUNDS ON CONDITION NUMBER: de (le) 9. 0137%6 | 


ALL VARIABLES IN THE MODEL ARE SIGNIFICANT AT THE 0.1000 LEVEL. 
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2. THE SECOND SAS MULTIPLE REGRESSION 


G@e7TONS LINESIZE=80;3 
DATA ONE; 
CMS FILEDEF DD1 DISK GDATA22 DATA A; 
mye ibe DDI; 
me UL KANG LOSS SPED TAMP; 
RS = RANGX*SPED; 
LS = LOSSxSPED; 
RL = RANGx*LOSS; 
PROC STEPWISE DATA=ONE; 
MODEL TAMP = RANG LOSS SPED RS LS RL / BACKWARD MAXR; 


SAS Gi 5o) Up oDAY ee UNE ws, 1989 
BACKWARD ELIMINATION PROCEDURE FOR DEPENDENT VARIABLE TAMP 
STEP 0 ALL VARIABLES ENTERED R SQUARE = 0.90984757 
CUP a= 7.00000000 
DE SUM OF SQUARES MEAN SQUARE id PROB>F 
REGRESSION 6 S6491.9759359S9 14415.32932267 94.20 Oy (HOM ela 
ERROR SE SomOwO to as a4 1 foo. OS7 05524 
HOTAL 62 95062.050469H4H0 
B VALUE Sip ee RACK LEVIS.” SCT pays F PROB>F 
INTCPT 109.07985496 
RANG =O. Ono 6 27 S002 S 505 17 Geresie@ 5027 0 iis 54 020013 
LOSS C2 les 7 164 C0029 772 3 4387.82533431 Zeros BES 80) (3yal 
Seed a OWA OMA, Leu ike Bl OO2 Fe aL Als 62Zee.co+lo262 heya is| 0.0001 
RS =O eOOCOUO1 33 Omer el evelem Hels 240.98633319 tee 0.2147 
LS 0.00000203 oooh eke hekekeh <x; ihe Fi sic) 70 Bekeh a ie) g tele! 0.0041 
RL = (2), 10,8) 8) OOO RSIS. 0.00000086 14723-42161660 A 2g Oieec 
BOUNDS ON CONDITION NUMBER: 7.239146, 206.9791 
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STEP 1 VARIABLE RL REMOVED RAR SQUARE = O2907960L6 


ChC2) a = 6.17240643 
DF SUM OF SQUARES MEAN SQUARE F PROB>F 
REGRESSION 5 86312.55431939 17262.51086388 112.46 0.0001 
ERROR 57 8749.49615001 153.49993246 
TORAL 62 95062.05046940 
B VALUE STD ERROR TYPE Tiss F PROB>F 
INTCPTF 112.32206019 
RANG =O O17 07 32 0.00264429 3008.87281840 19.60 0.000m 
LOSS 0.01136915 0.00168059 7024.89975456 u5.76 0.0001 
SPED -0.01705147 0.00269694 6136.0U920158 39.97 0.0001. 
RS -0.00000146 0200000105 292.39366190 17 cig 0.17249 
LS 0.00000206 0.00000068 1411.56803480 a, 20 0.003 Gm 
BOUNDS ON CONDITION NUMBER: 6.576496, 116.9728 
STEP 2 VARIABLE RS REMOVED R SQUARE = 0.90488434 
CUP) == 6.08301351 
DF SUM OF SQUARES MEAN SQUARE Bi PROB>F 
REGRESSION u 86020.16065748 21505.04016437 137395 0.0001 
ERROR 58 9041.88981191 155.89465193 
ISOM EP WE 62 95062.05046SH0 
B VALUE STD ERROR TYPE II SS F PROB>F 
INTCPT 118.23954426 
RANG =0-0107U358 O.001U7e57  155GGm72100 7605 q9,u3 0.0001 
LOSS 0.01154539 0.00168875 7286.435534u84 uG6.74 0.0001 
SPED =0.019970e76 0. COty Gm 1 OS ecco ee 125 86 0.0001 
ie 0.00000198 0.00000068 13290. 37593725 8.017 0.0051 
BOUNDS ON CONDITION NUMBER: u.72315, u5,35 


ALL VARIABLES IN THE MODEL ARE SIGNIFICANT AT THE 0.1000 LEVEL. 
SUMMARY OF BACKWARD ELIMINATION PROCEDURE FOR DEPENDENT VARIABLE 
VARIABLE NUMBER PARTIAL MODEL 


STEP RMVED IN Rxx2 Rxx2 GCE» F PROB>F 
1 RL 5 0.0019 0.9080 6.17241 1.1724 0.2838 
2 RS uy O20 si 0.9049 6.08301 1.9048 0.1729 
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SAS Pate Obra OE 13, 1989 

MAXIMUM R-SQUARE IMPROVEMENT FOR DEPENDENT VARIABLE TAMP 
meer 1 THROUGH 3 OMITTED 
STEP 4 VARIABLE LS ENTERED RB ace = 0.90488434 


OX = 5 SOs Osh) ail 
DF SUM OF SQUARES MEAN SQUARE F PROB>F 

REGRESSION > SOUC@G@eratea75S  A7iacs.04016437 $137.95 0.0001 
ERROR 58 9041.88981191 Poa.cg455193 
TOTAL 62 95062.05046940 

B VALUE STD ERROR PE VEIE, SE CROSS EF PROB>F 
myeCrT 118.23954426 
RANG -0.01474358 OOO / C7) sero ./ 1093605 99.43 0.0001 
LOSS Ol als ens ss) S CeCe lGs3s7 so 7286.435539484 46.74 0.0001 
See D POrmn01987067 Gog, 7 UCmmetoote 66957433 125.40 0.0001 
LS 0.00000198 0.00000068 eZee 99725 Sc47 OSG 0a1 
BOUNDS ON CONDITION NUMBER: ee ete Oi, EpeEO re) 


THE ABOVE MODEL IS THE BEST 4 VARIABLE MODEL FOUND. 


STEP 5 VARIABLE RS ENTERED R SQUARE = 0.90796016 
Gir) = 6.17240643 
DF SUM OF SQUARES MEAN SQUARE F PROB>F 
REGRESSION 5 86312.55431939 17262.51086388 112.46 hare veers 
ERROR 57 87 Homer o6 15001 153.4999324u6 
BOLPAL 62 95062.050469H0 
Boy AL oe STD ERROR BYPE iT SS a PROB>F 
MYFCPT 112.32206019 
RANG =O 0) 1 duge7.o 2 0.00264u429 S@@ene7e2sisu0 19.60 ale olantaral 
LOSS 0.01136915 0.00168059 7024.899754u56 45.76 O00 eM 
SPED -0.01705147 0.00269694 6136.04920158 39.97 0.0001 
RS -0.00000146 0.00000105 292.39366190 1.90 Oa 
LS 0.00000206 0.00000068 1411.56803480 a ,20 0.0036 
BOUNDS ON CONDITION NUMBER: 6.576496, dele .o72e 


THE ABOVE MODEL IS THE BEST 5 VARIABLE MODEL FOUND. 
weer 6 OMITTED 
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APPENDIX D. COMBINATORIAL OPTIMIZATION PROGRANI USING 


woe eu es 


x PROBLEM: 


ote 
«% 


t, 


Ses 


LAGRANGIAN RELAXATION 


eos. NTWRKS 


feslesleslestcsiesleslesics'es 


 eAiHORS: 


* SUBROUTINE: 


VARIABLES: 


A TACTICAL UNIT MUST MOVE FROM START POINT A, THROUGH 


AN INTERMEDIATE POINT B, TO OBJECTIVE © IN NO MORES iar 


TBAR TIME UNITS. EACH ARC/EDGE (U,V) IN THE NETWoR 


HAS A TRAVERSAL TIME TCCE) AND EXPOSURE/TACTICAL DIFFI- 
CULTY MOVEMENT ANTI=POTENTIAL OF TDCE). WE WISH TORR 


A PATH FROM A THROUGH B TO C WHICH HAS AN ACCEPTABLE 
TRAVERSAL TIME AND TACTICAL DIFFICULTY OR TO MINIMIZE 


TD SUBJECT TO TBAR. IT IS ASSUMED THAT THERE EXiS@ ae 


LEAST 1 FEASIBLE PATH FROM A THROUGH B TO C. 
RICHARD MILLER, MAJOR, USMC 
CHARLES H. SHAW, CAPTAIN, USA 
IAN KEITH, CAPTAIN, USMC 
DYKSTR - DIJKSTRA'S ALGORITHM 
UVCD - INPUT ARRAY OF EDGES (U,V),TIME COST TC(E), AND 
TACTICAL. DURE IC ULin vena. 
EDGES - NUMBER OF ARCS/EDGES (10000 MAXIMUM). 
X1,X2 - BOOLEAN VECTORS FROM X IN SUBROUTINE DIJKSTRA. 
OPTX - BEST VECTOR X1,X2,X3 FOUND THUS FAR. 
NODES - NUMBER OF NODES, I.E. , VERTICES (1100 MAXIMUM). 
S, T - START/STOP VERTICES FOR DIJKSTRA'S ALGORITHM. 
U, V - ALL ARBITRARY VERTICES. 
EPS - EPSILON, A SMALL NUMBER OR DIFFERENTIAL. 
TBAR - MAXIMUM TIME ALLOWED IN NETWORK (INPUT). 
PRED - PREDECESSOR OF A NODE. 
D - D(V) 'LABEL' ARRAY FOR DIJKSTRA'S ALGORITHM. 
DMIN - DMIN FOR DIJKSTRA'S ALGORITHM. 
EP - ENTRY POINT ARRAY. 


OUTDEG - QUTDEGREE OF NODE TO COMPUTE ENTRY POIND VAkRess 


ADJ - ADJACENT VERTICES ARRAY. 

LENGTH - EDGE LENGTH (TIME) ARRAY. 

TACDIF - EXPOSURE/TACTICAL DIFFICULTY ARRAY. 
TDMAX - LARGEST VALUE OF TACDIF(). 

OBJ - OBJECTIVE FUNCTION VALUE OF Pl IN DIJKSTRA. 
OBJMAX - OBJECTIVE FUNCTION MAXIMUM FOR P3. 

LAM - LAGRANGIAN MULTIPLIER IN SUBROUTINE DIJKSTRA. 
LAMMIN - TEMPORARY MINIMUM OF LAM. 

LAMMID - MIDPOINT OF LAMMIN AND LAMMAX. 

LAMMAX - TEMPORARY MAXIMUM OF LAM. 

ENDLAM - BEST VALUE OF LAM. 

TDLAMT - (TD + LAMBDA*T) ARRAY. 


SLOPE += SLOEE oe es : ee USN Bee NS pe Ses 


as Sas wo! teases e's mls we! tected eee! esis ste sles'esieste Pe sleslev 
VOTE SE SETS TEE TETE TESOL To ETE TOT TE TE TEC TE TE TS ETE TE TE TEAS TS ETE TE TOASTER TST TE TEC TE TEE COTE TE TS TET TDG TOG 


90 


Seseslesicsicstcsiesivsicsioe sleseslesicckatcleslescsleslesesccleslecescveslesesecesicsicslesesescvesteseccseslcvleveseseseseseste sk savesesk se 


gor T 


INTEGER UVCD(10000,4),X1( 10000) ,X2( 10000) ,OPTX( 10000) 
INTEGER EP( 10001) ,OUTDEG( 10001) 

INTEGER ADJ( 10001), LENGTH( 10001) , TACDIF( 10001) 
INTEGER EDGES,I,J,K,NODES,A,B,C,TDMAX,L,TIME,U,V 


REAL TBAR, LAMMIN, LAMMID, LAMMAX, OBJMAX, EPS, P1,ENDLAM 
REAL SLOPE1,SLOPE2,SLOPE3,0BJ1,0BJ2,0BJ3 


COMMON UVCD,EP,ADJ,LENGTH,TACDIF,EDGES ,NODES ,TBAR,A,B,C 
PARAMETER *CEPS = 0.0001) 
PEE ERS 


CALL EXCMS('FILEDEF 01 DISK PROJNET DATA B') 
CALL EXCMS('FILEDEF 02 DISK NETWORKS OUTPUT A’) 


* INITIALIZE ARRAYS TO ZERO 


5 


DO 5 1= 1, 10001 
OUTDEG( 1) = 0 


BEC) = 0 

ADJ(I) = 0 

LENGTH( 1) = 0 

TACDIF(1) = 0 
CONTINUE 


* READ IN THE DATA ONE ROW AT A TIME AND FIND 


“Bohs 


10 


Z0) 


NUMBER OF NODES AND MAXIMUM VALUE OF TD(E) 


READ(1,*) NODES,A,B,C,TBAR 
it 
TDMAX = 0 
READC1,*,END=20) (UVCD(1I,J), J=1,4) 
TDMAX = MAX(UVCD(1,4), TDMAX) 
= 
GO TO 10 
IF (NODES .GT. 1100) THEN 
STOP '2*** ERROR, NUMBER OF NODES EXCEEDS 1000 *7*' 
END IF 


* COMPUTE NUMBER OF EDGES IN THE NETWORK 


EDGES = I - 1 
IF (EDGES .GT. 10000) THEN 

STOP ‘*** NUMBER OF EDGES EXCEEDS SIZE DEFINED *?' 
END IF 


* COMPUTE OUTDEGREE FOR EACH NODE 


2 


DO=30 5) ei Deere 


OUTDEG(UVCD(I,1)) = OUTDEG(UVCD(I,1)) + 1 
30 CONTINUE 


* INITIALIZE ENTRY POINT ARRAY EPC) 


EPCI)e= = 

DO 40 I = 1, NODES 
EP(I+1) = OUTDEG(I) 
EPCI41 3 = -EPCUie = Leer 


40 CONTINUE 


%* COMPUTE THE HIERARCHICAL ADJACENCY LISTS FOR EACH SNODE EE. le 
* ADJ() AND LIST ALL THE NODE'S ADJACENT VERTICES AND CORRESPONDING 
* LENGTH OR TIME COST AND EXPOSURE/TACTICAL DIFFICULTY 


K=1 
DO 50 I = 1, NODES 
DO 51 J = 1, EDGES 
IF (UVCD(J,1) .EQ. I) THEN 
ADJ(K) = UVCD(J,2) 
LENGTH(K) = UVCD(J,3) 
TACDIF(K) = UVCD(J,4) 
K=K#1 
END IF 
IF (K .GT. 10000) THEN 
STOP ‘x ADJACENCY LISTS EXCEEDS SIZE DEFINED *%*' 
END IF 
51 CONTINUE 
50 CONTINUE 


PRINT *, ‘INITIALIZATION COMPLETE AND BEGINNING SEARCH' 


* FIND BOUNDS ON LAM 


LANMIN = 0.0 
LAMMAX = (NODES - 1) * TDMAX 
L= 0 


CALL DYKSTR( LAMMIN,X1,SLOPE1,0BJ1) 

WRITE(2,**) '‘LAMMIN=',LAMMIN,’ SLOPE1=' ,SLOPE1,' OBJ1=' ,OBJ1 
CALL DYKSTR( LAMMAX, OPTX , SLOPE3 ,OBJ3) 

OBJMAX = OBJ3 

WRITE(2,*) 'LAMMAX=',LAMMAX,' SLOPE3=' ,SLOPE3,' OBJ3=' ,OBJ3 
PRIAT *,' INITIAL DIJKSTRA COMPLETE ' 


70 LAMMID = (LAMMIN + LAMMAX)/2. 0 
IF (ABSC(CLAMMAX - LAMMIN) .LE. EPS) GOs7O 500 
Lb = Leet 


PRINT *; ITERATION {L, STARTING) 





80 


500 


85 


90 


25 


5518 
600 


CALL DYKSTR( LAMMID,X2,SLOPE2 ,OBJ2) 
WRITE(2,*) 'LAMMID=',LAMMID,' SLOPE2=',SLOPE2,' OBJ2=' ,OBJ2 


i CS hObE Za, Gi. §O}7 PHEN 
LAMMIN = LAMMID 
END IF 


ie SLOPEZ . LE. 0) THEN 
LAMMAX = LAMMID 
PEP Cend2 = Gl. ObUMAX) THEN 
OBJMAX = OBJ2 
ENDLAM = LAMNID 
DO 80 I] = 1, EDGES 
OPT SCL = x2 ( 1) 
CONTINUE 
END IF 
BND 
GO TO 70 


WRITE(2,*) 'WITH',L,‘ITERATIONS THE ' 
WRITE(2,85) OBJMAX 
FORMAT(' MAXIMUM OBJECTIVE FUNCTION VALUE FOR P2 IS ',F10. 3) 
iE = 0 
Pe Sa 
DO 90 I = 1, EDGES 
TIME = OPTX(I) * UVCD(I,3) + TIME 
P1 = OPTX(I) * UVCD(I,4) + Pl 
CONTINUE 
WRITE(2,*) 'THE SOLUTION FOR (TD*X) IS ‘',P1 
WRITE(2,*) 'EPSILON IS ',P1 - OBJMAX,', TIME ON THE PATH IS ‘',TIME 
WRITE(2,*) ‘AND THE LAST VALUE OF LAM IS ‘ ,ENDLAM 
WRITE(2,*) 'THE BEST PATH FROM NODE ',A,', THROUGH NODE ',B,’,. 
WRITE(2,*) ‘TO OBJECTIVE NODE ',C,'IS ALONG PATH: ' 
WRITE(2,95) A 
FORMAT( ‘NODE ',I5) 
U=A 
i 
ie Cue ne UVeD( I1.1)).AND. (OPTX(1)). BQ.91)) THEN 
WRITE(2,95) UVCD(I,2) 
ipCGeee. UVGD(1,2)) STOP THE END’ 
U = UVCD(I,2) 
GO TO 550 
ELSE 
=e eee 
GO TO 600 
END IF 


END 


wulenon'on'watontonton wnlon'eulon'enlon'en'eutonto s'¢ ate ate se se Js ale sle ole ale atenlcaleoate se s'e sc alent ontecatonton' enone se wleates aton'satenlontentecatententcuteatescatoatententen'enton'snton'e Pid when's 
¢ ex ever aver ex ev evaev ay et ev ev ex eh ev au eh Ai Geen Gh es en ew alt exe eh eh el au ese os eu ee eu ee a8 OF ae ah eV ee ee et es au at aN ah at el et aN aN ae eh aN aS aN ON eu ee ek 
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SUBROUTINE DYKSTR(LAM,X, SLOPE ,OBJ) 
INTEGER PRED(1100),1,J,K,S,T,U,V,VMIN, X( 10000) 
INTEGER UVCD(10000,4) ,EP( 10001) ,ADJ( 10001) , LENGTH( 10001) 
INTEGER TACDIF( 10001) ,EDGES ,NODES,A,B,C 
REAL D( 1100) ,DMIN, SLOPE, OBJ, LAM, TDLAMT( 10001), TINF,TBAR 
COMMON UVCD,EP,ADJ,LENGTH,TACDIF, EDGES ,NODES,TBAR,A,B,C 
PARAMETER (TINF = 1.0E15) 
*« SET TDLAMTC() FOR NEW LAM 
DO 300 I = 1, EDGES + 1 
TDLAMT(I) = TACDIF(I) + LAM * LENGTH(I) 
300 CONTINUE 
** RESET OPTIMAL PATH FOR NEXT SEGMENT 
DO 305 I = 1, EDGES 
MOD) = 
305 CONTINUE 
*« DIJKSTRA'S ALGORITHM FOR SHORTEST PATH FROM S TO ALL NODES 


A 
B 


5 
T 


Wool 


* INITIALIZE THE SHORTEST PATH AND PREDECESSOR ARRAYS 


400 DO 310) =e NODES 
PRED(I) = 0 
DCI) = TINF 
310 CONTINUE 


* INITIALIZE Sisk VERE es 


D(S) = 0.0 
PRED(S) = 8 


DO 320 I = 1, NODES 
DMIN = TINF 
DO 330 J = 1, NODES 
IF ((PRED(J).GE.0) . AND. (D(J).LT.DMIN)) THEN 
DMIN = D(J) 
VMIN = J 
END IF 
330 CONTINUE 
PRED(VMIN) = -PRED(VMIN) 
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DO 340 K = EP(VMIN), EP(VMIN+1)-1 
IF ((D(VMIN) + TDLAMT(K)) .LT. DCADJ(K))) THEN 
D(ADJ(CK)) = DCVMIN) + TDLAMT(K) 
PRED(ADJ(K)) = VMIN 
END IF 
340 CONTINUE 
320 CONTINUE 


mecONVERT PRED() TO PATH X() 


V=T 

360 I= 1 
Sere CC=PREDCV) . EQ. UVCD(1,1)) . AND. (V .EQ. UVCD(I,2))) THEN 

X(I) = 1 

V = -PRED(V) 

CLO ANON Sie 9, 
ELSE 

hai Oem REDS) sce TO 360 

=] oi 

COst0 276 
END IF 


Poel’ (5 .bQ. B) GO TO 390 


penow GET SHORTEST PATH FROM B TO C 


5 = B 
T= C 
GO TO 400 


SOMrUIE OBJECTIVE FUNCTION VALUE AND SLOPE FOR LAM 


390 OBJ = 0.0 
SLOPE = 0.0 
DO 450 I = 1, EDGES 
OBeie— (UVCDCT,4) + EAN UVCD(1,3))> = X(1) + OBJ 
SOR reee VGbCL 3) *~eXCL e+ oSLOPE 
450 CONTINUE 
2B = OBI S=shAMS TBAR 
SLOPE = SLOPE - TBAR 


END 


\O 
tay 


APPENDIX E. EXAMPLE MODEL OUTPUT 


The purpose of these example output plots is to demonstrate proper functioning of 
the model in terms of tactical movement decision logic and some sensitivity analvsis of 
the cognitive time scale values. Plots are appropriately labelled with all normal Tactical 
Optimum Path (TOP) output plots being expected value model runs unless otherwise 
noted. High risk averse plots use an input value of 30 minutes and extreme risk averse 
plots use an input value of 60 minutes rather than the expected value of 5.5 minutes. 

These verification model runs are based on the Lauterbach map sheet near Fulda in 
the central region of The Federal Republic of Germany. The base map 1s one version 
of the terrain relief map available in the CAMMS. The NATO vehicle represented is 
an MIAI Main Battle Tank moving as part of a deliberate attack against the specified 
Warsaw Pact systems shown in a deliberate defense. 

The sequence of numbers shown on the mapping represents the path points in the 
order specified by the user. The continuous line connecting these points 1s the resulting 
optimum path from the optimization. The other irregular polygon shapes represent the 
user input area of operations which appropriately restricts the network searched during 
the optimization. The set of optimum path points and large amounts of audit data are 
available in data files for use directly in combat models and analysis rather than output 


to a screen or map. 
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